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ABSTRACT 
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The  reaction  Be(d,n)  B  has  been  studied  by  tirae-of- flight 

techniques  at  bombarding  energies  of  3.5*  3.5*  and  5.5  MeV.  Contamin- 

l6  IT 

ant  neutron  peaks  due  to  the  0(d,n)  F  reaction  were  subtracted 
from  the  spectra  at  5.5  MeV  deuteron  energy.  Sixteen  levels  of  10B 
were  identified  up  to  7.56  MeV  excitation,  but  no  significant  evidence 
was  found  for  reported  states  at  5.18*  5.58,  6,Ho,  6,77  MeV.  The 
DWBA  theory  was  used  to  interpret  the  angular  distributions  for 
£-values  and  spectroscopic  factors,  and  plane  wave  theory  was  used  in 
cases  where  the  10B  levels  become  unbound  to  proton  emission.  The 
angular  distributions  for  the  states  at  3.59*  ^.77*  and  6,03  MeV  could 
not  be  satisfactorily  described  by  the  DWBA  theory.  The  spectroscopic 
factors  were  found  to  depend  most  sensitively  on  optical  model 
parameter  uncertainties,  and  also  on  whether  or  not  finite-range  and 
non-local  effects  were  included. 

Finally  in  an  appendix,  is  a  description  of  the  use  of  the  neutron 
time-of- flight  system  to  measure  target  thickness. 
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CHAPTER  I.  MOTIVATION  AND  THEORY 


1.1  Introduction 

This  thesis  reports  some  time-of -flight  studies  of  10B  states 
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excited  via  the  deuteron  stripping  reaction  Be(d,n)  Bf  at  bombarding 
energies  of  3*0,  3*5*  and  5*5  MeV.  The  purpose  behind  these  studies 
was  to  clarify  the  existence  of  several  weakly  excited  states  previ¬ 
ously  reported  for  ^B,  and  to  obtain  spectroscopic  information  for 
some  of  the  more  highly  excited  states,  notably  those  at  6.14,  6.57* 
6.88,  and  7.00  MeV. 

The  weak  states  of  interest  have  been  reported  at  excitations  of 
5*18,  5*58*  6.40  and  6.77  MeV.  We  shall  briefly  survey  the  existing 
evidence. 

The  5*18  level,  one  of  the  triplet  of  states  between  excitations 
of  5*10  MeV  and  5*20  MeV,  was  initially  interpreted  from  the  reactions 
6ttW}Y)10B  1),  2),  and  9Be(^)Y)10B  3).  With  a  centre-of- 

mass  width  estimated  at  (110  +  10)keV  ^),  observation  of  this  state  is 
difficult  although  it  has  been  confirmed  in  (d,n)  studies  at  3  MeV  ^). 

The  level  at  5*5 8  MeV  has  been  observed  in  (d,n)  work  by  Ajzen- 
berg  ^),  Hjalmar  ^),  and  in  a  (^Y)  experiment  by  Singh  ®),  but  the 
majority  of  evidence  is  negative  2*5*9°15).  jn  iat©r  high-resolution 
(j^Y)  studies  ^  ),  Y -transitions  associated  with  this  level  were  not 


detected. 


' 
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The  level  at  6*77  MeV  has  been  observed  in  (d,n)  work  only  ^),  and 
its  existence  is  largely  discredited  *  *  )* 

The  6*40  MeV  level  has  been  the  subject  of  some  controversy,  posi¬ 
tive  evidence  being  obtained  mainly  from  the  ^Be(d,n)^B  reaction 
6,7*9,15).  There  has  also  been  a  positive  result  from  the  ?Li(o^n)^B 
reaction  ^)*  On  the  other  hand  a  variety  of  reactions  leading  to 
states  have  yielded  no  trace  of  this  level* 


(i) 

9Be(3He,d)10B 

14.17) 

(iv) 

1XB(3He,(?010B 

(ii) 

10B(d,dB)10B* 

i2) 

11,13,19,20 

(iii) 

10B(p,p’  )10B* 

“> 

(v) 

2,21 

1 6 

It  must  be  pointed  out  that  0  target  contamination  can  give  rise 
to  spurious  lines  in  the  ^B  spectra*  For  example,  neutrons  from  the 
^0(d,n)^F  reaction  can  contaminate  forward  angle  neutron  peaks  corre¬ 
sponding  to  ^B  states  at  5»93  MeV  and  6*40  MeV.  The  state  at  0*5 
MeV  excitation  has  an  Ip  =  0  angular  distribution  of  neutrons  whose 
centre-of-mass  energies  approximately  correspond  to  those  from  the  ^®B 

state  at  6.40  MeV.  In  the  experiments  described  here,  the  contributions 
10  1  6  17 

to  the  B  spectra  from  the  0(d,n)  'F  reaction  were  estimated  by  sep¬ 
arate  measurements  using  an  oxygen-16  target. 

1*2  Theory 

The  plane-wave  (FWBA)  and  the  distorted-wave  (EWBA)  treatments  of 


the  stripping  reaction  mechanism  were  both  used  to  interpret  the  col¬ 
lected  data*  This  section  gives  a  description  of  the  theory  with  some 
discussion  of  the  validity  of  its  use  for  light  nuclei. 
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Stripping  theory  has  undergone  a  process  of  increasing  sophistica- 
tion  since  the  original  plane-wave  formulation  of  Butler  ).  Where 
the  deuteron  Instate  and  tensor  forces  are  ignored,  deuteron  stripping 
is  the  most  easily  described  of  such  rearrangement  collisions . 

In  the  case  of  (d,n)  stripping*,  Butler  proposed  that  a  measure¬ 
ment  of  the  neutron  angular  distribution  could  indicate  the  orbital 
angular  momentum  transferred  to  the  residual  nucleus.  This  leads  di¬ 
rectly  to  the  parity  change  in  the  reaction,  and  places  restrictions 

« 

on  the  spin  of  the  final  state .  For  a  zero- spin  target,  the  1-value 
measurement  specifies  the  spin  to  within  one  unit. 

With  a  reaction  A(d,n)B,  spin  conservation  requires  that 

4b  *  iA  +  Ip  + 1 p  *  £a  *  ip 

i.e.  | (min  ^  -^p  ^  ^a  ^  *^B 

where  s  spin  of  target  A, 

=  spin  of  residual  nucleus  B, 

1D  =  orbital  angular  momentum  transferred, 

■  intrinsic  spin  of  transferred  proton, 
jp  «  total  spin  of  transferred  proton. 

Parity  conservation  leads  to  the  equation 

nB  «  tt^-I)1? 

♦This  discussion  applies  equally  well  to  (d,p)  stripping  and  also  to 
the  inverse  pickup  reactions. 
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where  tt^  =  parity  of  the  ground  state  of  A* 

TTg  =  parity  of  the  final  state  of  B. 

Although  the  plane-wave  theory  enjoyed  some  initial  success  for 
light  target  nuclei,  it  failed  to  give  the  correct  absolute  cross- 
section,  often  overestimating  the  latter  by  an  order  of  magnitude  or 
more*  The  evidence  collected  from  stripping  on  medium  and  heavy  nuclei 
indicated  that  the  FWBA  was  not  wholly  successful  in  these  mass  regions, 
often  leading  to  ambiguous  /.-value  assignments.  Coupled  with  this 
fact  were  indications  that  the  simple  theory  did  not  give  the  correct 
dependence  of  the  differential  cross-sections  on  energy  and  Q-value. 

The  adaptation  of  the  theory  using  the  distorted  waves  Born  ap¬ 
proximation,  and  its  subsequent  application,  has  indicated  that  Coulomb 
and  nuclear  distortion  effects  are  always  important  in  stripping  re¬ 
actions.  The  plane-wave  treatment  is  never  a  good  approximation. 

The  LWBA  theory  allows  for  the  scattering  and  absorption  of  the 
incident  deuteron  before  stripping,  and  of  the  outgoing  nucleon  in  the 
exit  channel,  by  replacing  the  plane-waves  of  Butler* s  theory  by  elas¬ 
tic  scattering  waves.  The  latter  are  generated  by  Coulomb  and  optical 
model  potentials  which  reproduce  the  measured  elastic  scattering  from 
the  same  nucleus,  at  the  same  energy. 

There  have  been  many  excellent  reviews  of  the  modern  IWBA  theory 
23-26^  on]_y  a  brief  resume  of  the  version  used  for  the  present  cal¬ 
culations  will  be  given. 

The  differential  cross-section  for  A(d,n)B  stripping  can  be  written 
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*r.l  ^  kn  2JB  +  1  C2  2  S^B?2 

2  (2rrt>2)2  kd  2JA  +  1  i.m.j 


where  m^,  mn  are  the  reduced  masses  of  the  deuteron  and  the  neutron 
and  k^,  kn,  their  relative  wave  numbers,  C  is  the  Clebseh-Gordan 
isospin  coefficient 


^AtpMtA 


is  the  stripping  amplitude  ,  calculated  by  numerical  Integra** 
tion.  Explicitly 


0?  =  Jll 


-l 


Jzl  +1 


%D  (i&i'Xn)  (ip)  vnp/(5  0$d’£d)  dipd&i 


The  X(]£)  are  the  elastic  scattering  wave  functions  for  a  pair  of  part¬ 
icles  with  relative  momentum  k  and  separation  r.  The  plus  and  minus 


refer  to  their  asymptotic  behaviour,  is  the  wave  function  of  the 
proton  bound  to  the  target,  and  0^  is  the  deuteron  wave  function. 

The  residual  interaction  responsible  for  the  stripping  transition 
is  taken  as  Vnp,  the  neutron-proton  interaction.  More  accurately 


^residual  =  Vnp  +  ^nA  ~  Vng 

where  is  the  neutron-target  interaction  and  Vng  is  the  optical 
potential  in  the  final  channel.  One  usually  approximates  ~ 
although  this  is  not  strictly  true  since  is  real  and  VjjgLs  usually 
complex.  Additionally,  off-diagonal  elements  of  Vj^  which  correspond 
to  target  excitation  do  not  cancel  out,  and  one  hence  ignores  contri- 
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butions  from  the  two-step  processes  where  the  target  nucleus  is  excited 
by  inelastic  scattering  before  stripping  occurs 0 

Solutions  for  the  radial  parts  of^*,%“8  are  generated  numer¬ 
ically  from  the  Schrodinger  equation 


|"(r) 


2m 

t? 


[E 


E  +  V  (r)  +  V(r) 
c  o  m  o  c 


ti2  Id  f  d1  , 

2m  2 


0 


(1) 


where  V  (r)  is  the  Coulomb  potential  and  V(r)  the  optical  potential, 


m 


The  six-dimensional  integral  for  can  be  reduced  to  a  three¬ 


dimensional  form  by  using  the  zero-range  approximation  for  V  e  One  sets 

np 


Vnp  •  =  D*£(rnp) 

2  3/2 

where  D0  =  1023  x  10  MeV0fm  and  <f>^  is  the  Hulthen  wave  function a 

The  zero-range  approximation  essentially  means  that  the  deuteron  is 

absorbed  at  the  same  point  the  outgoing  nucleon  is  produced „ 

27 

The  "local-energy  "  approximation  )  allows  one  to  relax  the  zero- 
range  restriction  and  include  a  finite  range  for  the  n-p  interaction 0  In 

the  final  analysis,  this  approximation  includes  the  following  correction 

28 

factor  in  the  three-dimensional  integral  for  the  zero-range  case  )  s 


(  V  (r)  -  V  Jr)  -  V  (r. 
v  d  n  p 

d 

where  V^(r)  is  the  optical  potential  for  particle  i,  is  the  binding 
energy  of  the  deuteron,  and  B  =  7a,  with  $  “'1  *  0 „6  f m  . 


A(r)  =  (l  — 
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To  account  for  the  non-locality  of  nuclear  forces,  Perey^)  has 
proposed  that  local  wave  functions  (the  type  generated  by  equation  (1)) 
be  modified  as  follows % 


0NL  = 


1  +  “£  VL(r) 
.  2^2 


»l/2 


0! 


where  0^,  0^,  are  the  local  and  non-local  wave  functions  respectively, 
Vx,(r)  is  the  local  optical  potential.,  m  is  the  reduced  mass,  and  B  is 
the  range  of  the  non-locality^  usually  about  1  fm* 

Thus  non-local  and  finite-range  effects  are  included  by  multi¬ 
plying  the  radial  integrand  of  the  zero-range  case  by  four  form  factors* 
1*3  Discussion  of  the  Use  of  EWBA  Theory  for  Light  Nuclear  Reactions 
At  this  stage  one  can  say  that  stripping  reactions  on  medium  and 
heavy  nuclei  for  bombarding  energies  above  the  Coulomb  barrier  are 
reasonably  well  understood  in  terms  of  EWBA  theory* 

The  stripping  mechanism  for  light  nuclei  on  the  other  hand,  is  not 
as  clearly  understood*  In  fact  one  can  advance  several  reasons  why  one 
would  not  expect  standard  stripping  theory  to  be  very  successful  for 
light  nuclei  especially  at  low  bombarding  energies* 

Some  difficulties  stem  from  compound  nucleus  effects  which  are 
important  in  light  nuclei  because  of  their  lower  level,  density*  At  low 
energies  the  compound  nucleus  can  usually  exist  in  a  series  of  sharp 

states  well  separated  in  energy*  If  the  energy  spread  in  the  incident 

spacing 

beam  is  less  than  the  mean  level/  then  th®  scattering  cross-section  will 


show  characteristic  variations  unaccounted  for  by  the  optical  model* 
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Compound  elastic  scattering  is  indistinguishable  from  the  direct 

shape  elastic  scattering  described  by  the  optical  modelo  Hence  optical 

% 

parameters  for  low  energy  scattering  in  a  light  nuclear  region  may  be 
in  error.  Also  for  very  light  nuclei ,  exchange  reactions  may  take 
place  sufficiently  often  to  render  the  optical  model  invalid  except  for 
small  scattering  angles. 

Inelastic  neutron  scattering  will  also  contribute  in  cases  where 
(d,n)  stripping  is  being  studied.  It  has  been  noted  however  ^)#  that 
despite  strong  fluctuations  in  differential  cross-sections  near  a  com¬ 
pound  resonanceB  the  angular  shapes  vary  little.  This  may  indicate 
that  the  competing  compound  amplitude  for  many  reactions  in  light  nu¬ 
clei  is  small.  In  the  latter  studies  of  the  ^CCd^nJ^N  reaction  for 
energies  1.50  MeV  to  MeV  it  was  found  that  the  shape  of  the  angu¬ 

lar  distributions  averaged  over  bombarding  energy  is  probably  given  by 
the  direct  amplitude  alone. 

Statistical  model  (Hauser-Feshbach)  calculations  performed  by 
Siemssen  et  al  for  the  ^Be(d[,n)^^B  reaction  at  energies  between 
1*5  MeV  and  3*0  MeV  were  not  wholly  successful.  They  concluded  that 

1 

the  statistical  model  calculations  either  highly  overestimated  the 
compound  amplitude  ,  or  that  th©  statistical  model  was  not  valid  for 
the  open  channels  considered.  In  th©  light  of  these  statements  it 
was  felt  that  no  significant  information  could  be  gained  by  performing 
statistical  model  calculations  for  the  ^Be(dt(n)‘L^B  reaction  at  5*5  MeV, 
A  recent  study  of  the  latter  reaction  at  ?.00  MeV  bombarding  energy 
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where  the  appropriate  optical  parameters  were  known,,  indicated  that  the 
EWBA  theory  was  moderately  successful  in  reproducing  the  main  stripping 
peaks,  at  least  for  transitions  to  the  first  four  states  of  ^B# 

One  of  the  difficulties  of  using  DWBA  theory  for  a  reaction  such 
as  ^Be(d,n)^B  is  that  the  extreme  j-j  coupling  model  does  not  hold  for 
very  light  nuclei.  Generally  above  mass  16,  j-j  coup3.ed  wave  functions 
are  used  in  the  stripping  formulation#  However  it  is  known  that 
the  spin-coupling  scheme  changes  from  the  L-S  extreme  at  the  beginning 
of  the  Ip- shell  to  the  j-j  form  as  the  shell  fills#  Unfortunately,  for 
most  of  the  studied  transitions  more  than  one  j-value  is  allowed  from 
angular  momentum  conservation#  Hence  there  is  a  chance  that  a  mixture 
of  different  j-values  will  be  present.  In  the  calculations  reported 
here,  the  extreme  j-j  coupling  model  was  employed,  and  a  single  strip¬ 
ping  amplitude  allowed  for  in  any  one  calculation  0 
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CHAPTER  II  o  EXPERIMENTAL  DETAILS 


2 . 1  Introduction 

In  fast  neutron  spectroscopy,,  the  most  precise  measurement  of 
neutron  energy  is  by  means  of  a  i±me-af- flight-  spectrometer  which 
determines  the  energies  of  neutrons  by  individually  timing  their  flight 
over  an  accurately  measured  patho  Sine©  the  amplitude  of  a  signal  from 
a  neutron  scintillation  detector  is  rarely  related  to  the  particle 
energy  in  an  unambiguous  way,  a  pulse  height  distribution  alone,  is 
insufficient  for  neutron  energy  determination  * 

The  basic  function  a  time-of-flight  spectrometer  performs  is  to 
accurately  fix  two  time  points? 

(i)  the  zero  or  start  time  of  the  neutron  from  the  target 
(ii)  the  time  the  neutron  is  detected  at  the  end  of  the  flight  path# 
From  this  time  measurement,  the  neutron  velocity  is  known  and  hence 
the  energy,, 

These  related  time  problems  are  best  solved  by  pulsed  beam  tech¬ 
nique  which  in  recent  years  has  developed  into  the  most  powerful 
method  available  for  fast  neutron  spectroscopy  in  the  nano-second 
time  scale. 

Beam  pulses  of  short  duration  in  comparison  to  the  flight  times 
to  be  measured,  strike  the  target  and  initiate  the  reaction.  The 
time  zero  signal  marking  the  arrival  of  the  ion  burst  at  the  target 
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is  obtained  by  passing  the  particle  beam  through  a  hollow  cylindrical 
capacitor  before  the  target.  The  small  charge  pulse  induced  is  ampli¬ 
fied  and  appropriately  shaped. 

The  corresponding  stop  signal  is  obtained  from  the  neutron  de¬ 
tector  which  usually  consists  of  an  organic  scintillator  optically 
coupled  to  a  photomultiplier  tube.  The  time  between  beam  pulses  is 
sufficiently  long  that  the  lowest  energy  neutron  of  interest  can  be 
detected  before  the  arrival  of  the  next  pulse. 

The  start-stop  pulses  are  then  fed  t©  a  fast  time-to-amplitude 
converter  and  the  corresponding  output  signal  analysed  by  a  multi¬ 
channel  pulse  height  analyser. 

The  attractive  features  of  the  pulsed  beam  technique  include  high 
resolution,  high  detection  efficiency  and  applicability  to  a  wide  range 
of  neutron  energies.  The  use  of  high-speed  tunnel  diode  switching 

circuits  has  introduced  considerable  improvement  in  spectrometer  reso- 

| 

lution  in  the  past  few  years.  Improvements  in  beam  pulsing  techniques 
have  also  contributed. 

In  the  next  two  sections,  a  description  of  the  main  experimental 
system  is  given.  Section  2.2  deals  with  the  essential  features  of  the 
beam  pulsing  equipment  installed  in  this  laboratory,  and  section  2.3 
describes  the  electronics  associated  with  the  time-of-flight  spectro- 

i* 

meter. 

2.2  Beam  Pulsing  Equipment 

The  accelerator*  and  associated  beam  pulsing  equipment  delivered 

*6.00  MeV  Van  de  Graaff  Accelerator,  Model  CN,  High  Voltage  Engineering 
Corporation,  Burlington,  Mas sachus setts,  U.S.A. 
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to  the  target  assembly  high  intensity  deuteron  bursts  of  optimum  pulse 
length  0c5  ns  at  a  repetition  rate  of  1  MHa, 

Several  beam  pulsing  methods  have  been  developed  in  the  pasts 

(i)  chopping  the  continuous  ion  beam  from  a  Van  de  Graaff  by  means 

30 

of  a  beam  deflector  and  slit  )? 

(ii)  pulsing  the  ion  beam  before  acceleration 

(iii)  compression  of  the  pulsed  beam  produced  in  (ii)  after  accelera- 
32 

tion  J  )• 

For  deuterons  of  several  MeV  energy  the  most  successful  method  is 
a  combination  of  pre-acceleration  pulsing  and  post-acceleration  com¬ 
pression  of  this  pulsed  beam  using  a  Mobley  bunching  magneto  High  peak 
currents  (^ma)  can  be  obtained  with  very  short  pulse  durations  (<  Ins) 
and  background  radiation  is  much  reduced  over  method  (i)» 

Such  a  system  has  been  installed  at  the  University  of  Alberta 
Nuclear  Research  Centre  and  a  simplified  schematic  of  the  equipment  is 
given  in  figo  lo  A  detailed  treatment  of  this  system  with  special  em¬ 
phasis  on  the  theory  of  the  Mobley  magnet  buneher  can  be  found  in  the 
PhoDo  thesis  of  W,Go  Davies*^) <>  A  brief  description  only  will  be 
presented  here# 

At  the  high  voltage  Van  de  Graaff  terminal 9  the  R,Fo  ion  source 
delivers  a  d9c«  current  of  several  railliamps  of  deuterons*  The  beam 
is  focussed  on  a  l/l6  inch  aperture  and  swept  across  the  aperture  in 
an  elliptical  pattern  by  a  pair  of  R,F,  deflection  plates  set  at  right 
angles  to  each  other  in  space 9  with  the  RoF„  voltages  90°  out  of  phase. 
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FIG  I  *  ESSENTIAL  FEATURES  OF  BEAM  PULSING  EQUIPMENT 


' 


14 


This  pulsing  kit  produces  a  triangular  pulse  with  a  full-width  at  half- 
maximum  (FWHM)  of  between  10ns  and  15ns  at  a  repetition  rate  of  1MHz. 

The  deuterons  ar©  accelerated  and  focussed  at  the  accelerator  base 
then  refocussed  through  a  90°  mass-energy  analysing  magnet*  Energy 
stabilisation  of  the  accelerator  is  achieved  by  a  corona  servo¬ 
mechanism  controlled  by  signals  from  slits  at  the  focus  of  the  analysing 
magnet* 

A  magnetic  quadripole  lens  refocuses  the  beam  at  a  point  between 
the  plates  of  the  R*F*  deflection  system  for  the  Mobley  magnet  buncher. 
Phase  information  ©f  the  ion  puls©  is  supplied  to  the  10  MHz  R.F.  oscil¬ 
lator,  from  a  hollow  cylindrical  electrode  which  surrounds  the  beam. 

A  90°  Mobley  magnet  then  delivers  a  time- compressed,  space-focussed 
pulse  to  the  target  assembly* 

The  basic  idea  behind  a  Mobley  bunching  system  is  to  apply  a  trans¬ 
verse  momentum  modulation  to  the  initial  10ns  ion  pulse  and  then  to 
sweep  this  deflected  pulse  through  a  suitable  uniform  field  magnet  in 
such  a  way  that  all  particles  in  the  pulse  arrive  at  the  target  at 
almost  the  same  time  in  a  well-focussed  spot*  The  beam  bunching  factor 
for  the  present  installation  is  of  the  order  of  12  using  a  90°  magnet 
with  a  radius  of  75”  •  An  optimum  pulse  duration  of  0*5ns  is  possible. 

The  main  factor  which  limits  the  bunching  performance  in  terms  of 

33 

time  dispersion,  is  a  property  called  the  beam  emittance  •  The 
latter  gives  a  measure  of  the  transverse  dimensions  and  the  angular 
divergence  of  the  pre-bunched  puls©,  and  is  a  constant  for  motion  in  a 
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uniform  magnetic  fieldo  Beam  quality  and  the  focussing  properties  of 
the  quadrupole  lens  are  the  most  significant  factors  governing  the 
beam  emittance. 

2.3  Time-of-Flight  Electronics 

During  the  period  in  which  part  ©f  this  experimental  work  was 
carried  out,  a  fast  digital  computer  and  memory  display  unit*  became 
available,  allowing  large  quantities  of  data  t©  be  processed  within  a 
reasonable  time.  Apart  from  the  computer  control  of  data  collection, 
the  main  time-of-f light  electronics  used  in  the  primary  studies  of  the 
^Be(d,n)^^B  reaction,  were  similar  to  the  system  operating  under  pro¬ 
gram  control.  Consequently  the  latter  system  will  be  described  in 
this  section. 

The  SDS  920  computer  has  8192  24-bit  words  of  memory  of  which  4096 
are  available  for  data  storage.  A  cathode  ray  tube  display  unit  oper¬ 
ating  under  program  control  allows  a  two-  or  three-dimensional  display 
of  data  collected  in  memory. 

Connected  on-line  i.e.  as  a  stored  program  multidimensional  ana¬ 
lyser  3^»35)  ^  .^e  computer  allowed  two  ma  jor  improvements  to  be  made 
in  the  existing  time-ef -flight  system. 

1)  The  versatile  general  purpose  kicksorter  program  (GPKS),  written  by 
Dr.  W.K.  Dawson  could  perform  a  real-time  correction  for  time-shift  ) 
in  the  detection  time  of  the  fast  signals  from  the  neutron  detector.  A 


♦Scientific  Data  Systems  920,  1649  Seventeenth  Street,  Santa  Monica, 
California. 
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FIG. 2,  BLOCK  DIAGRAM  OF  THE  TIME~OF«FLIGHT  EL ECTRONICS 
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description  of  the  time-shift  correction  is  given  in  Appendix  A® 

2)  Immediately  a  parti eular  run  was  completed,  a  large-scale 
visual  record  ©f  the  collected  spectra  could  be  obtained  via  the 
c®r®t®  unito  This  enabled  electronic  drifts  or  poor  beam  quality  to 
be  quickly  diagnosed® 

A  schematic  block-diagram  of  the  electronics  is  given  by  fig®  2, 
Essentially  there  are  two  parallel  time-of-f light  systems®  The  first 
records  the  neutron  time  spectra,  the  second  provides  a  monitor  of 
the  first  record,  allowing  corrections  to  be  made  for  a  fluctuating 
count  rate,  due  to  changing  background  or  beam  conditions® 

37 

Considerable  experimental  analysis  of  the  on-line  system  )  has 
indicated  that  at  low  neutron  energies,  the  best  time-resolution  is 
obtained  when  the  leading  edge  of  the  fast  signal  from  the  principal 
neutron  detector  is  used  t©  specify  time-aero®  The  stop  signal  is 
provided  by  the  conventional  cross-over  detection  of  the  delayed  beam 
pick-off  pulse® 

Four  signals  are  taken  from  the  principal  neutron  detector  which 
consists  of  a  liquid  scintillator*  coupled  to  a  Phillips  XP1040 
photomultiplier  tube® 

The  fast  time  signals  are  taken  from  the  I4th  dynode,  inverted, 
and  fed  to  a  start  discriminator  shaper  designed  by  Dr®  GoC®  Neilson. 
The  leading  edge  of  this  signal  provides  the  time-zero  for  the  time- 

*Ne213  in  a  cylindrical  quartz  container,  3rol/2"  diameter  by  3/4*'  deep® 
Nuclear  Enterprises,  Ltd®,  Winnipeg,  Manitoba.® 
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(jQ 

to-amplitude  converter  (TAG)  which  is  based  on  a  design  by  DoL.  Wieber  ) 
The  stop  signals  for  both  time-of-flight  systems  are  obtained  from 
the  cylindrical  beam  sensor  located  30  cm  before  the  target „  This 
start-stop  logic  ensures  that  the  stop  discriminator  will  trigger  only 
when  a  neutron  is  detected 0 

Signals  from  the  anode  and  dynode  13  of  the  Phillips  tube,  provide 
the  imput  to  a)f-n  discriminator  of  the  pulse-overlap  variety  0  Sig¬ 
nals  from  dynode  12  have  a  three-fold  purposes 

I*  to  provide  a  pulse  height  spectrum  for  the  time-shift  correction; 

2<>  to  provide  a  gate  signal  for  the  time-to-amplitude  converter; 

3«  to  provide  a  coincidence  with  the  neutron  signals  from  the  Y~ 
eliminator. 

The  pulse  height  spectrum  is  fed  to  a  128  channel  analogue-to- 
digital  converter  (ADC  b)  and  the  output  from  the  TAC  to  the  1024  chan¬ 
nel  converter  ADC  A» 

With  the  time-shift  correction,  a  system  dead-time  of  the  order  of 
700  y^s  allows  GPKS  to  process  a  maximum  of  150Q  events  per  sec.  During 
the  time  GPKS  is  storing  an  event,  input  to  the  computer  is  inhibited 
by  the  coincidence  interrupt  unit.  The  latter  demands  a  triple  coinci¬ 
dence  between  the  following  three  signals  before  input  to  the  SDS  920 
occurs  s 

1,  an  address  reset  pulse  from  ADC  A; 

2,  an  address  reset  pulse  from  ADC  B; 

3*  a  logic  pulse  from  coincidence  unit  C, 
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Signals  1,2  ensure  that  both  ADC“s  accept  amplitude  pulses  from 
the  experiment  only  when  both  are  free  to  do  so#  and  signal  3  requires 
that  the  accepted  pulses  represent  a  true  coincident  pair.  A  logic 
pulse  from  the  interrupt  unit  to  the  computer  then  indicates  that  the 
next  two  pulses  on  the  data  lines  are  to  be  treated  as  a  coincident 
pair*  A  detailed  description  of  this  interrupt  unit  is  contained  in 
an  internal  report  by  Do A.  Gedcke^  )* 

The  monitor  detector  consists  of  a  1"  Naton  phosphor*  mounted  on 
an  RCA  8575  photomultiplier  tube*  The  fast  signals  are  taken  from  the 
anode,  inverted,  and  detected  in  a  fast  tunnel  diode  cross-over  pick- 
off  circuit*  The  input  to  the  TAC  is  gated  by  suitably  shaped  signals 
from  the  13th  dynode*  Low-energy  neutrons  and  ]f-ray  pulses  are  elim¬ 
inated  by  a  single- channel  analyser* 

The  output  from  the  monitor  TAC  is  analysed  and  stored  by  a  256 
channel  pulse-height  analyser,  the  contents  of  which  can  be  later 
dumped  into  the  computer  memory  under  control  of  GPKS* 

2*4  The  9Be(d,n)1QB  Reaction 

Using  the  preceding  time-of- flight  system,  this  reaction  was 
studied  at  three  deuteron  energies  3»00®  3°50o  5»54  MeV*  The  average 
current  of  deuterons  on  the  target  was  between  2,0^ a  and  2,5yu.a*  The 
9 Be  target  consisted  of  a  thin  film  of  Beryllium  metal  vacuum  deposited 
on  a  l/32  inch  thick  Molybdenum  foil*  Various  flight  paths  between 
4*0m  and  6.3m  were  used,  but  the  majority  of  data  was  collected  at  a 
♦Nash  and  Thompson  Ltd.,  Hookrise  South,  Tolworth,  Surrey,  England. 
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flight  path  of  6.3m  and  a  deuteron  energy  of  5*5^  HeV.  For  this  deu- 
teron  energy  data  were  collected  at  5°  intervals  from  0°  to  35°,  and  for 
15°  intervals  from  45°  to  150°  (laboratory  angles). 

The  experiments  at  3*0  MeV,  3*5  MeV  deuteron  energy  constituted  a 
further  search  for  the  broad  5*18  MeV  level  of  ^B.  For  these  energies 
measurements  were  made  at  the  following  laboratory  angles:  0°,  15°, 

30° ,  45°,  60°. 

2.5  The  ^0(d.n)~^F  Reaction 

With  the  purpose  of  estimating  the  contribution  of  this  reaction 
to  the  ^B(d,n)^B  spectra  at  5*54  MeV,  a  separate  study  of  this  reaction 
was  carried  out  using  an  oxygen-16  target.  The  latter  consisted  of  a 
thin  film  of  tantalum  pentoxide  (Ta20^),  deposited  on  a  tantalum  foil. 

As  far  as  possible,  experimental  conditions  were  the  same  as  those  used 
for  the  ^Be(d,n)^B  reaction  at  5*5^  MeV,  and  measurements  were  made  at 
the  same  laboratory  angles.  As  a  check  against  the  similarity  of  the 
two  experiments,  data  for  the  ^Be(d,n)^B  reaction  were  collected  at 
the  same  time,  for  laboratory  angles  of  0°  to  45°  inclusive. 
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CHAPTER  III .  DATA  ANALYSIS 


3.1  Introduction 

The  analysis  and  interpretation  of  the  collected  data  can  be 
divided  into  two  parts.  The  first  part  contained  in  Chapter  III  deals 
with  the  reduction  of  the  raw  ^Be(d,n)^B  spectra  to  yield  excita¬ 
tions  and  relative  angular  distributions  of  neutrons.  The  second  part 
described  in  Chapter  IV,  presents  evidence  for  the  four  weak  states 
under  study  and  describes  the  interpretation  of  the  neutron  distribu¬ 
tions  for £ -values  and  spectroscopic  factors  using  plane-wave,  and 
distorted  wave  theories  of  stripping, 

A  typical  ^Be(d,n)^B  spectrum  for  5*5^  MeV  deuteron  energy  is 
indicated  in  fig,  3*  The  range  of  neutron  energy  co  vered  by  this  spec¬ 
trum  is  from  about  1,5  MeV  at  the  neutron  cut-off,  to  10  MeV,  corre¬ 
sponding  to  the  ground  state  of  ^B,  For  these  spectra  an  overall  time 
resolution  of  1,3  ns  was  achieved,  corresponding  to  an  energy  resolution 
of  85  keV  for  6  MeV  neutrons  over  a  6,3  m  flight  path.  At  this  energy, 
sixteen  neutron  groups  were  identified  corresponding  to  ^B  excitations 
up  to  7*56  MeV, 

The  data  at  3*0  MeV,  3*5  MeV  were  collected  after  the  neutron  and 
Y-ray  background  had  been  reduced  with  additional  shielding.  A  much 
improved  time  resolution  of  0,85  ns  was  achieved  corresponding  to  an 
energy  resolution  of  about  50  keV  for  6  MeV  neutrons  over  a  6.3  flight 
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path*  A  3.0  MeV  spectrum  for  a  laboratory  angle  of  15°  is  indicated 
in  fig.  4. 

3.2  Extraction  of  Excitations  and  Relative  Angular  Distributions 

Due  to  the  considerable  quantities  of  data  collected  in  present 
day  nuclear  physics  experiments,  reduction  to  meaningful  information 
requires  the  use  of  fast  computers,  if  the  analysis  is  to  be  completed 
in  a  reasonable  time. 

A  flow  diagram  for  the  data  analysis  of  these  experiments  is  pre¬ 
sented  in  fig.  5.  In  all } eight  separate  computer  programs  were 
employed,  the  majority  of  which  were  written  in  this  laboratory. 

The  in  itial  task  was  identification  of  known  neutron  peaks  in  the 
spectra  using  a  relativistic  kinematics  program  written  by  T.B.  Grandy. 
From  the  known  excitations  and  bombarding  energies  the  program  supplied 
the  corresponding  neutron  flight  times  for  a  6.3  m  flight  path,  thus 
enabling  a  time  scale  to  be  constructed  for  each  spectrum.  The  times 
for  the  unknown  peaks  could  then  be  found. 

The  program  involved  with  the  estimation  of  the  mean  excita¬ 
tions  and  the  Gaussian  analysis  of  the  spectral  peaks  were  written  by 

33 

W.G.  Davies,  and  are  fully  reported  in  his  Ph.D.  thesis  ). 

Firstly,  a  general  purpose  least-squares  program  was  employed  to 
fit  a  linear  or  quadratic  time  calibration  line  for  all  the  spectra, 
using  the  times  for  the  known  peaks.  This  line  was  computed  at  all 
angles  because  at  some  angles  slight  differences  were  introduced  by 


electronic  drifts. 
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Fig. 5.  Flow  diagram  of  data  analysis. 
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Using  the  coefficients  of  these  lines,  and  the  flight  times  for  a 
given  neutron  peak  at  all  angles,  best  mass  and  Q-values  were  estimated 
again  using  the  least  squares  method.  On  deciding  that  a  given  neutron 
peak  did  correspond  to  the  ^Be(d,n)^B  reaction,  the  correct  mass 
was  inserted  and  a  mean  ^B  excitation  calculated. 

In  cases  where  neutron  peaks  were  not  resolved,  or  contaminated 
from  other  reactions,  this  estimation  was  postponed  until  after  the 
Gaussian  analysis  of  the  spectra. 

The  areas  of  the  neutron  peaks  were  calculated  using  a  weighted 
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non-linear  Gaussian  analysis  program.  Using  the  least-squares  method  ), 
this  program  fits  a  function  of  the  following  form  to  a  set  of  N  peaks: 


N 

y  =  background  +  2 

i  =  1 


Aj_  exp 


“(x  -  xi) 


2 


Here  y  denotes  the  count  in  channel  x.  The  background  function  could 
have  a  variety  of  forms  —  linear,  quadratic,  or  exponential, -but  in 
practically  all  cases,  a  quadratic  form  was  employed.  For  each 
Gaussian  curve,  the  three  parameters  A±t  6\t  and  Xi  denoting  respec¬ 
tively  the  peak  height,  the  standard  deviation,  and  the  mean,  had  to 
be  estimated.  Up  to  100  channels  of  data  with  ten  free  parameters 
could  be  dealt  with  in  any  one  calculation. 

In  the  determination  of  errors  for  the  above  quantities  Ai» 

,  x^,  the  program  includes  both  the  statistical  error  and  the 
chi-squared  fitting  error  ^). 

Explicitly,  the  mean  square  error  of  an  estimated  parameter  ^  \ 
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(e.g.  mean  ,  standard  deviation  etc.)  from  its  true  value 
is  given  by 

<K-  =  jgL  .<W 

where  (i)%^  =  I4  W±  (jj±  -  B  i)^  with  the  weight,  yU.-^  the  measured 
value  on  the  Gaussian  curve,  and£j_  the  theoretically  estimated  value; 

(ii)  N  is  the  number  of  measurements  and  p  is  the  number  of  free 
parameters ; 

(iii)(0(<4)f  is  the  mean-square  deviation  in  o estimated  from  the 
least-squares  analysis. 

A  comparison  between  the  mean  square  error  in  and  its  mean- 
square  deviation  indicates  the  importance  of  non-statistical  errors. 

If  systematic  errors  can  be  ignored,  a  value  of  0(?/(N-p)  significantly 
greater  than  unity  would  imply  that  the  Gaussian  function  used  in  the 
analysis  is  wrong. 

The  assumption  of  a  Gaussian  form  for  the  spectral  peak  is  open 
to  some  question,  but  in  most  cases  this  has  proved  to  be  a  reasonable 
assumption.  Peaks  with  large  areas  tended  to  be  slightly  asymmetrical, 
but  the  areas  calculated  using  the  program  were  always  within  %  of 
the  values  calculated  by  other  means,  e.g.  summation  by  hand  with 
subtraction  of  a  constant  background.  The  danger  of  using  the  auto¬ 
matic  method  is  that  it  may  consistently  under-  or  over-estimate  the 

area  thus  introducing  systematic  errors.  A  better  method  has  been 
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described  recently  by  Tepel  ). 
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The  Gaussian  fit  for  all  peaks  of  interest  at  one  angle  supplied 

the  fitting  information  (means,  standard  deviations  etc.)  which  was 

then  used  to  calculate  the  necessary  fitting  information  for  the  other 

fifteen  angles.  With  this  control  information,  the  estimation  of  peak 

areas  etc.  could  be  placed  on  a  semi-automatic  basis. 

In  order  to  calculate  angular  distributions  using  these  peak 

areas,  the  efficiency  of  the  principal  neutron  detector  as  a  function 

of  neutron  energy  had  to  be  determined.  This  was  done  by  studying  the 

reaction  ^Li(p,n)^Be  for  proton  energies  between  3*0  MeV  and  5*0  MeV. 

At  0°,  this  proton  energy  range  gave  rise  to  neutrons  of  energies 

between  1*31  MeV  and  3*33  MeV.  The  neutrons  were  detected  both  by  the 

principal  neutron  detector  and  a  " flat-response”  long  counter  of  the 

41 

McKibben  type  )•  The  detectors  were  placed  side-by-side  at  the  longest 

flight  path  available  (6.3m),  in  such  a  position  that  they  subtended 

the  same  angle  with  respect  to  the  proton  beam  direction. 

The  long  counter  detects  neutrons  corresponding  to  the  ground  and 

7 

first  excited  states  of  Be.  The  excited  state  counts  were  removed 
from  the  total,  by  dividing  by  the  ratio  (1  +  I^/Iq)  measured  by 
Bevington  et  al  ) .  The  total  ground  state  counts  detected  by  the 
principal  neutron  detector  were  then  normalised  to  the  corrected  McKibben 
counts  for  each  proton  energy,  and  the  resulting  efficiency  curve  is 
indicated  in  fig.  6.  The  solid  curve  is  a  least-squares  fit  to  the 


data  using  the  function 

Efficiency  =  A 
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The  symbols  A,  E0,  n  are  parameters  of  the  fit,  and  CT np  is  the  semi- 


2X2  \  -1 


<3"np  =  3tt  ( 1.206  En  +  (-1.86  +  0,09415  En  +  0.0001307  En  ) 


After  a  study  of  the  monitor  spectra  taken  at  each  angle  it  was 
decided  to  normalise  the  ^Be(d,n)^^B  spectral  areas  to  the  total  counts 


in  the  first  200  monitor  channels. 

For  those  peaks  free  of  contaminants,  relative  angular  distribu¬ 
tions  were  then  calculated  using  a  program  written  by  T.B.  Grandy. 

This  program  normalised  the  areas  to  the  monitor  sums  and  made  correc¬ 
tions  for  the  relative  efficiency  of  the  neutron  detector  and  for 
neutron  attenuation  in  the  gold  target  backing  (l/l6n  thick).  The 
angular  distribution  program  combines  errors  from  the  following  sources: 

1.  statistical  errors  from  fitting  process 

2.  errors  estimated  for  efficiency 

3«  errors  in  monitor  sums. 

In  a  manner  exactly  similar  to  the  above  description,  relative 

l6  17 

angular  distributions  were  obtained  from  the  0(d,n)  'F  data,  for 
neutron  groups  corresponding  to  the  ground  and  first  excited  states 
of  ^F.  The  two  angular  distributions  are  shown  in  fig.  7. 

Then,  essentially  by  inverting  the  calculations  in  the  angular 
distribution  analysis,  the  actual  neutron  contributions  to  the 
^Be(d,n)10B  spectra  from  the  1^0(d,n)1^F  reaction  were  estimated. 
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FIG. 6  THE  RELATIVE  EFFICIENCY  OF  THE  NEUTRON  DETECTOR. THE 
SMOOTH  CURVE  IS  A  LEAST-SQUARES  FIT  TO  THE  DATA  USING  THE 
INDICATED  FORMULA  DUE  TO  BAME^. 
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FigeT*  The  relative  angular  distributions  of  the  0(d9no;  F9 
^OCdjn^)  reactions  at  Ed»5e5^M*V«  The  relative  yields  are 
given  in  Appendix  B  on  page  73 * 
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a  n  q 

These  estimates  were  made  both  for  the  initial  7Be(d#n)  B  data  at  5*54 

"L6  17 

MeV,  and  for  the  data  collected  during  the  study  of  the  0(d,n)  'F 
reaction. 

In  this  fashion  it  was  possible  to  construct  a  complete  angular 
distribution  for  the  5*93  MeV  state  of  ^Bt  and  to  estimate  upper  limits 
for  excitation  of  the  controversial  6.4 0  MeV  state. 


b  -v  •  tot  a  r  -  •;  S»c  e:  >dfT 


. 


. .  1 .  J  :  . 

. 


■  '  > 


. 

. 


. 


. 


33 


CHAPTER  IV.  DATA  INTERPRETATION 

4.1  Results  of  Excitation  Analysis  for  Weak  States 

A  careful  study  of  all  the  measured  spectra  has  revealed  no  sig¬ 
nificant  evidence  for  the  broad  level  at  5*18  MeV.  The  spectrum  in 
fig.  4  represents  a  run  with  the  best  time  and  energy  resolution  which 
has  been  achieved  to  date  with  the  spectrometer  used  in  this  laboratory. 
Although  the  most  modern  evidence  appears  to  confirm  the  existence  of 

the  5*18  MeV  level,  no  conclusive  result  can  be  presented  here. 

"L6  17 

The  effects  of  removing  the  0(d,n)  'F  contamination  from  the 

^Be(d,n)^B  spectra  indicate  that  the  peak  occurring  at  6.40  MeV  is 

17 

probably  due  entirely  to  excitation  of  the  0.5  MeV  level  of  F.  The 
largest  uncertainties  of  the  contamination  were  estimated  at  not 
greater  than  10$  for  forward  angles. 

This  conclusion  is  at  variance  with  the  recent  results  of  Bussino 
and  Smith  who  report  an £p  =  o  angular  distribution  for  a  "^B  state 
at  excitation  (6.35  ±  0.05)  MeV.  Their  DWBA  fit  to  this  neutron  dis¬ 
tribution  was  for  angles  between  0°  and  35°  and  no  account  was  taken 
1  17 

of  possible  0(d,n)  'F  contamination. 

Negative  results  are  also  reported  here  for  the  states  at  5*58  MeV 
and  6.77  MeV.  The  upper  limits  of  possible  contributions  due  to  the 
four  weak  states  are  presented  in  Table  1. 
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The  upper  limits  of  states  in  at  excitations  of  5*18» 
5»58>  6.40,  6*77  MeV.  The  5«l8  MeV  limit  is  relative  to  the 
yield  for  the  (5.11  +  5*16)  MeV  doublet.  The  5*58,  6.40, 

6.77  MeV  limits  are  relative  to  the  yield  for  the  ^®B  ground 
state.  Both  the  areas  of  (5.H  +  5*16)  MeV  and  0.0  MeV 
states  are  set  to  1.00. 


Ed 

3.50  MeV 

5.54  MeV 

State 

5.18 

5.58 

6.40 

6.77 

Angle 

0° 

<  0.03 

<  0.03 

<  0.04 

<  .003 

1 5° 

<  .02 

20° 

<  .02 

<  0.04 

<  .005 

25° 

<  .03 

o 

o 

<  .03 

<  .02 

<  .03 

<  .009 
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<  .02 
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4.2  Angular  Distribution  Interpretation 

The  angular  distributions  of  fourteen  neutron  groups  were  resolved 

for  excitations  up  to  7.00  MeV,  and  were  interpreted  for ^-values  and 

spectroscopic  factors  using  both  the  plane-wave  stripping  theory  ****), 

45 

and  the  IWBA  code  of  Smith  )•  The  latter  makes  provision  for  approx- 

28 

imating  the  finite-range  n-p  interaction  ),  the  non-locality  of  nu- 

29  46 

clear  forces  7),  and  also  the  deuteron  polarizability  ). 

The  optical  potentials  for  the  deuteron  and  neutron  channels  took 

the  following  forms? 

U(r)  =  - 

f(r)  =  ( 

and  either  g(r)  =  f(r)  (Volume  absorption) 


Vf(r)  +  iWg(r)j  where 


1  +  exp  ^ 


R' 
sC 


-i 


with  R  =  R0A' 


1/3 


or  g(r)  =  4a  4^  (derivative-type  surface  absorption). 

dr 

The  derivative-type  form-factor  is  to  be  preferred  however,  since  it  is 

known  that  the  most  important  contributions  to  stripping  reactions  come 

64 

from  the  nuclear  surface  ). 

The  interior  wave  function  of  the  captured  proton,  E,  ,  is  calcu- 

J/'P 

lated  by  starting  with  an  asymptotic  solution  and  then  working  inward 
numerically  toward  the  origin.  The  potential  well  has  a  real  Woods- 
Saxon  part  -Vf(r),  no  imaginary  part,  but  does  include  a  derivative- 
type  spin-orbit  well  of  depth  6  MeV.  The  real  part  V  is  changed  by 
trial-and-error,  until  the  radial  solution  at  the  origin  has  less  than 


1-  •'  F?j  K  j  O  £  t(V"  4  '  is  *Mlo 

".ol  ^  f  j  .>1  eiJ.t 


%  „ 

-  ii.  J  »»  bj  w  ti  otcr  ■  bo  i*.-»  ,o  :Ii  >  nr  Jfcnu&t  rv&v/  loiteJ  i 

' 


36 


l/lOOOth  of  the  value  for  the  solution  at  r  s  Rp  +  4ap. 

The  solutions  for  both  the  deuteron  and  neutron  radial  wave  func¬ 
tions  F^„  s  are  estimated  by  starting  with  zero  at  the  origin^lO"-^ 
for  r  =  0*15  fta»  and  then  working  outward  numerically  using  equation 
(1)  (Section  1.2) ,  until  the  calculated  values  match  with  prescribed 
asymptotic  forms. 

Calculations  were  performed  both  with  and  without  a  cut-off,  using 
the  zero-range  and  the  finite-range  approximations  with  non-local 
effects  included* 

Optical  parameters  for  the  elastic  scattering  of  7  MeV  deuterons 
from  ^ Be  ^ )  f  and  14  MeV  neutrons  from  ^B  ^),  led  to  EWBA  differential 
cross-sections  in  poor  agreement  with  the  5° 5  MeV  data.  A  subsequent 
parameter  search  was  undertaken  for  EWBA  fits  to  the  ^B  ground  state 
angular  distribution  using  the  following  options? 

1*  A  zero-range  calculation  with  Woods-Saxon  form  factors  for  the  real 
and  imaginary  optical  well  depths  (ZR). 

2.  A  zero-range  calculation  with  Woods-Saxon  form  factors  for  the  real 
parts  of  the  optical  potentials  and  with  derivative-type  form  factors 
for  the  imaginary  parts*  These  calculations  were  repeated  with  radial 
cut-offs  of  2*5  fm„  3»0  3<>5  fin.  (ZRCO)* 

3*  A  finite-range  calculation  with  form  factors  as  in  2  but  no  cut-off  (FR). 
4*  A  finite-range  calculation  with  form  factors  as  in  2  but  with  a 
cut-off  of  3«°  fm  (FRCO). 

5*  A  finite-range  calculation  with  form  factors  as  in  2,  no  cut-off. 
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but  including  non-local  effects  using  a  range  of  p  =  O.85  fm  for  the 
non-locality  in  the  three  channels  (FRNL). 

The  purpose  behind  the  us@  ©f  these  options  was  to  investigate 
their  effects  on  the  magnitude  and  shape  of  the  DWBA  differential  cross- 
sections  calculated  using  the  same  optical  parameters  for  all  the 
options  * 

Initially  only  option  1  was  available  for  calculations  and  an  ex¬ 
tensive  parameter  survey  yielded  many  parameter  combinations  which  gave 
a  good  overall  fit  to  the  ground  state  angular  distributions.  For  this 
survey,  the  radii  and  diffuse ness  parameters  were  fixed  at  the  following 
values ; 

R0d  =  1.80  fm  s  Rop  =  1.5  fra  8  Ron  =  1.3  fm 

ad  =  0.7  fm  =  ap  g  an  s  0.5  fm. 

This  large  value  of  the  deuteron  radius  has  been  used  in  previous 

49 

calculations  at  3  MeV  ).  The  following  search  proceedure  was  adopted: 

(a)  Vn  was  varied  alone  between  40  MeV  and  70  MeV  in  5  MeV  steps; 

(b)  was  varied  alone  between  45  MeV  and  120  MeV  in  5  MeV  steps; 

(c)  Wn,  W^  were  initially  held  fixed  at  5  MeV,  10  MeV  respectively; 

(d)  a  given  fit  was  then  improved  by  adjustment  of  Vn»  Vd*  Wn,  Wd  in 
1  MeV  steps. 

A  list  of  favourable  sets  of  Vn,  Vd,  Wn,  Wd,  is  indicated  in  table  2. 

The  program  employed  an  externally  controlled  increment  proceedure  to 
determine  optimum  parameter  combinations .  The  theoretical  and  experi¬ 
mental  distributions  were  superposed  in  a  logarithmic  plot  and  a  fit 
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TABLE  20  Potential  sets  yielding  acceptable  fits  to  the  B(d,n0)  B 


distribution  using  a  zero-range  DWBA  calculation * 


NOo 

-V 

n 

-w 

n 

-vd 

-wd 

(  MeV  )  -Vp  =  l+5MeV 

1 

kk 

6 

50 

9 

R  =  lfl30fm  ;  a  =  0  0  5fm 

2 

b$ 

k 

kQ 

8 

on  n 

3 

55 

5 

79 

9 

R  ,  «  l«80fms  a,  =  0,7fm 

k 

55 

k 

77 

2 

oa  d 

5 

57 

k 

80 

8 

R  =  1  o  50fm;  a  =  0  07fni 

6 

57 

k 

79 

10 

°P  P 

7 

57 

5 

78 

11 

for  all  sets  listed 

8 

60 

5 

81 

8 

9 

60 

k 

80 

10 

opposite  c 

10 

60 

k 

79 

11 

11 

62 

k 

82 

8 

12 

62 

5 

81 

10 

13 

— 

6k 

k 

82 

9 

TABLE  3o  Potential  sets  yielding  acceptable  fits  to  the  ^Be(d,no)10B 


distribution  using  a  finite-range  DWBA  calculation  with 
non-local  effects  included 0 


o 

o 

£3 

-V 

n 

-W 

n 

-vd 

-W  (  MeV  ) 

n 

-Vp  =  l+0MeV 

Ik 

55 

9o5 

50 

606 

R 

on 

=  lo35fm  i  a  =  O055fm 
n 

15 

50 

k,0 

55 

10o0 

R  A 

od 

■  l0^8fm  i  a,  =  0063fm 

d 

16 

55 

6e0 

62 

10o0 

R 

op 

«  l«U8fm  g  a  =  0c63fm 

P 
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all  sets  listed  opposite* 
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was  tested  visually 0  A  value  of  chi-squared  was  calculated  for  each 
pair  of  distributions  and  this  also  gave  a  measure  of  the  goodness- 
of-fit. 

When  options  2,  3»  4,  5  became  available  with  the  Smith  code,  a 
parameter  search  was  undertaken  using  option  5°  For  this  case  the 
radii  and  diffuseness  parameters  were  fixed  at  the  following  values? 

Rod  =  1.48  fra  s  R0p  =  1.48  fm  i  Ron  =  I.35  fra  47 ) 
ad  =  O.63  fm  s  ap  =  O.63  fm  s  a„  =  0.55  fm 
Initially  the  values  of  Wn,  were  held  fixed  at  9*5  MeV  and 
6.6  MeV  respectively*  (These  values  were  also  taken  from  refs*  4?,  48). 

The  search  proceedure  of  option  1  Was  then  carried  out.  With  this 
survey,  fewer  parameter  sets  were  found  which  gave  the  quality  of  the 
fits  obtained  with  option  1.  None  of  the  sets  of  table  2  gave  a  good 
overall  fit  for  the  ground  state  distribution  although  the  main  strip¬ 
ping  peak  was  generally  reproduced*  Some  of  best-fit  sets  determined 
from  this  survey  are  listed  in  table  3» 

Using  the  potential  sets  of  tables  2  and  3  with  the  appropriate 
options,  I1WBA  cross-sections  were  generated  for  all  the  resolved  neu¬ 
tron  distributions  corresponding  to  bound  levels.  States  in 
with  Q- values  of  less  than  -2*226  MeV  became  unbound  to  proton  emission, 
and  the  DWBA  code  used  was  unable  to  describe  stripping  to  virtual 
levels.  (There  is  a  DWBA  theory  for  such  stripping  but  it  has  been 
used  in  only  one  published  case  to  date  ^)).  As  a  point  of  interest, 
computational  difficulties  arose  with  the  Smith  code  for  Q- values  less 
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than  -2.00  MeV.  This  was  found  to  be  due  to  the  numerical  range  allowed 
by  the  IBM  7040  computer  which  performed  the  calculations » 

The  angular  distributions  which  could  not  be  described,  corresponded 
to  10B  excitations  at  6*57  MeV  (Q  =  =2.21  MeV),  6.88  MeV  (Q  =  -2.52  MeV), 
and  7.00  MeV  (Q  =  =2.64  MeV). 

Satchler  and  Smith  *^)  have  both  advised  that  the  Q- values  for 
these  states  be  adjusted  until  a  bound  condition  is  obtained  (Q  =-2.00 
MeV).  One  can  argue  heuristically  that  provided  the  "captured11  proton 
remains  significantly  below  the  Coulomb  plus  centrifugal  barrier,  its 
wave  function  looks  presumably  much  like  a  bound  state  save  for  a  slow 
decay.  Other  than  seeming  reasonable,  there  is  no  physical  justification 
of  this  proceedure. 

Regarding  the  EWBA  fits,  it  was  found  that  a  parameter  set  yield¬ 
ing  a  good  overall  fit  for  the  ground  state  distribution,  generally  did 
not  accurately  describe  the  excited  state  distributions  away  from  the 
main  stripping  peak.  Single  parameter  sets  however  resulted  in  unam¬ 
biguous  ^-values  for  the  following  states  (MeV)s 

0.0,  0.72,  1.74,  2.15,  5.11,  5.16,  5. 93,  6.14. 

These  ^-values  were  confirmed  by  plane-wave  calculations  performed  con¬ 
currently  using  a  " stripping"  radius  of  5  fin  for  all  the  calculations. 

The  results  of  EWBA  calculations  for  ten  bound  states  are  presented 
in  figs.  8,  9,  and  10,  using  the  single  parameter  set  listed  in  fig.  8. 
Superposed  for  each  state  are  the  curves  resulting  from  options  2,  3,  5* 
The  plane-wave  calculations  for  the  states  at  6.57,  6.88,  7.00  MeV  are 
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0,0,  0.72,  1.74,  2,15,  AND  3.59  MeV  FOR  E  ■  5.54  MeV.  THE  SMOOTH  CURVES 
ARE  DWBA  CALCULATIONS  WITH  THE  PARAMETERS  LISTED  f  USING  THE  OPTIONS 

DESCRIBED  IN  THE  TEXT. 


:  ■  ’  •  .8.m 


k2 


.CENTRE  OF  MASS  ANGLE(degrees) 

FIG. 9.  ANGULAR  DISTRIBUTIONS  OF  THE  NEUTRONS  LEADING  TO  THE  10B  STATES  AT 
h, 77,  5.11,  5-16  MeV'FORE  ■  5-51*  MeV.  THE  SMOOTH  CURVES  ARE  DWBA  CALCULATIONS 
USING  THE  PARAMETERS  AND  OPTIONS  LISTED  ON  FIG. 8, 


■ 


CENTRE  OF  MASS  ANGLE  (degree  s)  10 

FIG.  10.  ANGULAR  DISTRIBUTIONS  OF  THE  NEUTRONS  LEADING  TO  THE  •  B  STATES  AT 

i 

5.93f  6.03,  6,14  MeV,  FOR  Ed=  5.54  MeV.  THE  SMOOTH  CURVES  ARE  DWBA  CALCULATIONS 
USING  THE  PARAMETERS  AND  OPTIONS  LISTED  ON  FIG. 8. 


EXPERIMENT  e  =  5.5U  MeV 


kk 


CENTRE  OF  MASS  AN  GLE  ( degrees) 

FIG.  11.  ANGULAR  DISTRIBUTIONS  OF  THE  NEUTRONS  LEADING  TO  THE  10B  STATES  AT  6,57,  6.88,  AND  7.01  MEV. 
THE  SMOOTH  CURVES  ARE  PLANE-WAVE  (BUTLER)  CALCULATIONS  USING  A  ’STRIPPING ’  RADIUS  OF  5.0  FM  THROUGHOUT. 
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shown  in  fig*  XI*  The  results  of  EWBA  calculations  for  these  states 
proved  unsuccessful* 

The  *^B  states  at  excitations  3«59v  4*77  and  6.03  MeV  had  anoma¬ 
lous  neutron  distributions  which  defied  the  direct  reaction  interpreta¬ 
tion  (fig  So  8-10) 

15 

The  3°59  MeV  level  has  been  previously  assigned  J?p  =  1  ),  but 

for  this  Si  -value  the  LWBA  code  generated  a  peak  at  the  wrong  angle* 

The  4*77  MeV  level  has  always  given  trouble  in  stripping  reactions 
14*15).  xt  has  recently  been  assigned  Jn  -  3+  both  experimentally  ^)9 
and  theoretically  from  effective  interaction  calculations  in  the  lp 
shell  -^)*  However  neither  Xp  ”  1*  nor  ^p  =  3  (values  required  from 
parity  conservation)*  was  successful  in  generating  the  measured  shape, 
which  appears  devoid  of  the  familiar  direct  reaction  structure  (fig.  8)* 
In  fact  as  2  was  the  most  successful  lvalue*  although  this  also  was 
unsatisfactory*  If*  as  Kurath  has  suggested  ^)*  this  state  is  due  to 
a  two-particle  excitation  of  the  nucleus*  then  the  single-particle 
description  employed  here  is  inadequate* 

The  6.03  MeV  level  has  a  spin  of  4+  ^)*  a  value  which  would  re¬ 
quire  £,p  =  3  as  a  minimum.  No  /.-value  between  0  and  3  resulted  in 
satisfactory  EWRA  calculations*  Further  analysis  of  ^Be(d*n)^B  data 
at  3*5  MeV  (collected  July  17*  1964  in  this  laboratory)  did  not  elim¬ 
inate  this  anomaly. 

The  best  plane-wave  estimates  for  the  states  at  6*57*  6.88*  7<>00 
MeV  gave  JLp  =  3*  fp  s  1,  (p  =  1  respectively*  However  the  last  two 
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TABLE  4.  MEASURED  STATES  OF  10B 


This  Experiment 

Previous  Work* 

Ex  (MeV  ±  keV)  £p 

ip 

Ex  (meV  £  keV) 

j" 

T 

0.0 

1 

1 

0.0 

3+ 

0 

0.72  ±  10 

1 

1 

0.7173 

±  0.8 

1+ 

0 

1.74  ±  10 

1 

1 

1.740 

±2 

0+ 

1 

2.15  ±  10 

1 

1 

2.154 

±  3 

1+ 

0 

3-58  ±  10 

(1) 

1 

3.585 

±  4 

2+ 

0 

4.77  ±  10 

(2.3) 

>2  15) 

4.774 

±3 

3+  35) 

0 

5.11  ±  12 

0 

0 

5.114 

±  ^ 

(2") 

0 

5.17  ±  14 

1 

1 

5.166 

±  * 

2+ 

1 

5. 183 

±8 

1+ 

0 

5*58 

6,8) 

5o93  ±  10 

1 

1 

5«923 

2+ 

0 

6.03  ±  12 

6.029 

±4 

4+ 

6.14  ±  10 

(2) 

(2) 

6.133 

±  ^ 

"  t  r 

(0) 

6.35 

±  50  15) 

6.57  ±  10 

(3) 

1 

6.566 

±  6  , 

6.77 

6) 

6.89  ±  15 

(i) 

6.884 

1" 

0 

7.00  ±  12 

(i) 

7.00 

7.431 

±10 

2” 

0 

(7.468 

±  10) 

(2+) 

7.48  ±  15 

7.479 

±2 

(2‘) 

7.56  ±  25 

- 

7.561 

±1 

0+ 

1 

68 

♦All  information  taken  from  ref.  )  unless  otherwise  marked. 
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-values  must  be  adopted  with  caution  because  of  probable  contamination 
of  the  neutron  distributions  from  ^C(d,n;j_)^N  reaction* 

The  Rvalues  from  this  set  of  experiments  ,  together  with  the  mean 
■^B  excitations  and  errors ,  estimated  with  the  least-squares  code,  are 
con^ared  with  previous  measurements  in  table  4, 

The  fact  that  a  ground  state  best-fit  set  of  parameters  did  not 
yield  a  good  overall  fit  for  the  excited  state  distributions,  prompted 
an  investigation  of  the  effects  ©f  varying  the  values  of  Vn  and  Wn  for 
the  excited  states.  Generally  a  given  fit  could  be  improved  by  adjust¬ 
ing  Vn  and  Wn  alone.  However,  there  did  not  appear  to  be  any  system¬ 
atic  correlation  of  the  neutron  parameters  with  excitation  energy, 

49 

Siemssen  et  aL  )  have  suggested  that  the  optical  parameters  for  such 
a  light  nucleus  might  depend  ©n  the  configuration  of  the  outer  nucleons 
which  presumably  varies  from  state  to  state. 

However,  considering  the  uncertainties  of  compound  nucleus  contri¬ 
butions,  and  the  inadequacies  of  th©  stripping  formulation  for  this 
mass  number  and  bombarding  energy,  it  is  hardly  surprising  that  one  set 
of  parameters  does  not  completely  describe  the  excited  state  distri¬ 
butions  , 

On  the  basis  of  calculations  performed  on  the  relative  angular 
distribution  alone,  it  is  difficult  to  compare  th©  merits  of  options 
1  to  5*  For  a  given  potential  set  with  the  same  radial  form  factors, 
only  small  changes  in  the  calculated  angular  shapes  are  detectable. 
Presumably  option  5#  with  finite-range  and  non-local  effects  included 
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represents  the  most  accurate  calculation* 

The  use  of  a  radial  cut-off,  advocated  by  Satchler  to  elimi¬ 
nate  erroneous  contributions  from  the  nuclear  interior  is  a  most 
drastic  solution  although  it  appears  to  improve  a  zero-range  fit  for 
back  angles. 

From  an  extensive  series  of  calculations  performed  on  the 
^Ca(d,p)^Ca  reaction  ^),  it  is  known  that  the  inclusion  of  finite- 
range  and  non-local  effects  tends  to  dampen  contributions  from  the 
nuclear  interior,  thus  reducing  the  need  for  a  cut-off.  However,  the 
non-local  correction  tends  to  enhance  the  tail  of  the  bound  state  wave 
function,  so  that  although  interior  contributions  are  diminished,  the 
cross-section  magnitude  is  unchanged  from  calculations  performed  with¬ 
out  a  cut-off. 

As  an  illustration  of  these  effects,  the  different  deuteron  radial 
wave  functions  for  the  ground  state  transition,  calculated  using  options 
3,  5  are  indicated  in  fig.  12. 

In  their  study  of  the  ^Ca(d,p)^xCa  reaction  from  7.00  MeV  to 

57 

12.00  MeV  deuteron  energy,  Lee  at  al  )  came  to  the  conclusion  that 
the  best  agreement  with  experiment  could  be  achieved  by  using  a  deu¬ 
teron  potential  of  real  depth  120  MeV,  roughly  twice  the  nucleon 
potential.  Recently,  Greider  and  Strobel  have  advocated  using  a 
deuteron  potential  equal  to  the  sum  of  the  neutron  and  proton  poten¬ 
tials.  For  nuclei  where  rnp  ^  ,  they  indicate  how  such 

a  deuteron  potential  more  accurately  justified  the  use  of  Vnp  as  the 
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residual  interaction  in  stripping. 

Although  Greider  and  Strobe!  cs  conditions  are  not  met  with  in  the 
9Be(d,n)^B  reaction,,  some  finite-range  EWBA  calculations  were  performed 
for  deuteron  potentials  between  80  M©V  to  140  MeV  ^V^)  A  tolerable  fit 
to  the  ground  state  distribution  was  obtained  using  the  FRNL  option  for 
a  deuteron  potential  of  about  110  MeV.  However  this  aspect  of  the 
calculations  was  not  pursued  in  any  detail  and  no  calculations  for  the 
excited  states  were  made. 

4.3  Calculation  of  Spectroscopic  Factors 

In  the  case  where  only  one  lvalue  and  j-value  are  important  in  a 
stripping  transition,  the  spectroscopic  factor  can  be  defined  by  the 
normali s ation 

(Cr(6))ex(j=  $i‘  ( 

However  for  only  three  of  the  measured  states— 0o0,  1.74,  5«H  MeV— 
were  the  j- values  known  explicitly  (j  =  3 /2#  3/2  o  l/2  respectively). 

For  the  others ,  j  s  {  +  l/2,  or  j  s  /-  l/2,  were  allowed.  Since  the 
bound  proton  has  a  spin-orbit  well  of  depth  6.00  MeV,  the  ambiguity  in 
j  will  be  reflected  in  uncertainty  in  the  absolute  EWBA  cross-sections, 
and  hence  in  the  spectroscopic  factors.  For  the  states  where  uncer¬ 
tainty  in  j  existed,  calculations  were  performed  with  both  j-values , 
and  this  led  to  EWBA  cross-sections  for  a  given  state  differing  by  a 
maxi mum  of  7$.  All  the  spectroscopic  factors  reported  here  have 
j  =  H  +  l/2.  Since  relative  cross-sections  only  were  measured  the 
spectroscopic  factors  have  been  normalised  to  a  ground  state  value  of 
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TABLE  5.  RELATIVE  SPECTROSCOPIC  FACTORS 


♦Potential! 
Set  ! 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

State 

Gnd 

1 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0  ] 

L.O 

1.0 

1.0 

1.0 

0.72 

1 

2.76 

2.79 

2.79 

2.77 

2.78 

2.75 

2.75 

2.73 

2.72  2 

i.75 

2.69 

2.70 

2.65 

1.74 

1 

1.83 

1.79 

1.74 

1.82 

1.82 

1.91 

1.97 

2.04 

2.13  2 

>.18 

1.95 

1.84 

1.99 

2.15 

1 

.78 

.71 

.66 

.69 

.64 

.63 

.63 

.60 

.58 

.61 

.59 

.59 

.58 

3.59 

1 

.43 

.37 

.6? 

.64 

.72 

.67 

.66 

.70 

.65 

.60 

.67 

.63 

.64 

4.77 

1 

.56 

.48 

.55 

jih. 

0  TT 

.62 

.58 

.56 

.65 

.58 

.54 

.45 

.48 

.45 

4.77 

2 

.45 

.40 

.29 

.36 

.64 

.79 

*72 

1.08 

1.06  1.07 

.84 

.88 

.80 

5.H 

0 

.59 

.54 

.51 

.58 

.54 

.61 

.66 

.46 

.53 

.58 

.41 

.46 

.45 

5.16 

1 

.55 

o62 

.74 

.75 

082 

.99 

1.00 

lol6 

1.18  1.17 

1.12 

1.12 

1.0  7 

5.93 

1 

.81 

.73 

.69 

.63 

.74 

.79 

.81 

1.10 

1.13  1.17 

1.18 

1.05 

1.17 

♦♦6.14 

2 

1 .83 

.26 

.36 

.32 

.43 

.4? 

.44 

.60 

.65 

.61 

.53 

.69 

.63 

♦The  number  indicates  the  potential  set  in  Table  2. 
♦♦A  spin  of  3  has  been  assumed  for  this  level* 
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TABLE  60  Relative  spectroscopic  factors  for  all  ^B  states  up  to 
601^  MeV  excitation9  using  the  potential  sets  listed  in 
TABLE  3o 


Potential  No. 

Option 

lb 

ZRCO 

lb 

FRCO 

lb 

FR 

14 

FRNL 

15 

FRNL 

16 

FRNL 

STATE 

ip 

s£ 

Gnd0 

1 

loO 

loO 

loO 

1,0 

1,0 

1,0 

0cT2 

1 

2c9B 

2o96 

2o93 

2o65 

2  068 

2,72 

loT^ 

1 

10^5 

iM 

lo5^ 

l0bh 

lo6l 

1.U8 

2015 

1 

0c^6 

0  046 

0o51 

0ch9 

0,56 

0,53 

3o59 

1 

0o2^ 

0o2U 

0028 

0,27 

0o33 

0029 

k.77 

1 

OcOT 

0o0T 

0e09 

0C08 

Ooll 

0o09 

U.77 

2 

Oc36 

0o38 

0cb2 

0M 

0o33 

0.39 

5  oil 

0 

0o39 

0  c  50 

0,70 

o  062 

O063 

0,55 

5d6 

1 

0>9 

0  c  49 

0o6l 

0c59 

O065 

0,5^ 

5*93 

1 

Oo59 

Qo59 

0oj8 

Oc75 

0,73 

O067 

*6clk 

2 

Oo50 

Oc  50 

0  oT2 

O063 

0o58 

0,55 

*  A  spin  of  3  has  been  assumed  for  this  state  when 
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TABLE  T«  Comparison  of  measured  relative  spectroscopic  factors  with 
previous  experimental  and  theoretical  values,, 


Spectroscopic  Factors 


1 

STATE 

£p 

Plane  wave 

values 

Average  DWBA 

values 

Previous  DWBA 

values1^) 

Theoretical 

S(lp)2B  S(lp)Pot* 

Gnd 

i 

1,0 

1,0 

1,0 

1,0 

1.0 

0,72 

i 

2,75 

2,77t ,1C 

2,30 

1,38 

1,7^ 

i 

1.27 

l,78i,2C 

1,00 

1,96 

2,15 

i 

0,6l 

o,6o±«ii 

0,1+1 

0,U2 

3,59 

i 

0,29 

0,51-«20 

0,1+5 

0,09 

0.32 

**.77 

i 

0,06 

0,02 

0,01 

U.77 

2 

o.i+o 

0,68±,30 

0.1*3 

5,11 

0 

1,66 

0.5U±,13 

1.30 

5,16 

1 

0.1*7 

0,82^25 

0,16 

0.23 

5.93 

1 

0,66 

0,81+^15 

0,51+ 

0.35 

6,  ih 

2 

0,30 

0,521,22 

*  The  spectroscopic  factor  for  a  lp  nucleon  is  the  number  of  lp  nucleons 

in  10B  times  the  squares  of  the  coefficients  of  fractional  parentage 

10  ^/2  1/2 
for  the  B  states  in  quest ion 9 summed  over  the  lp  and  lp 

components o  Kurath  has  done  this  for  two  lp~shell  residual  interactions 

-  the  {8-l6)Pot  and  the  (6~l6)2B  described  in  ref,  ^),  For  both 

interactions,  the  ground  state  of  10B  has  an  absolute  theoretical 

spectroscopic  factor  of  1,2, 
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1.0. 

Relative  factors  were  calculated  for  all  the  potential  sets  of 
tables  2  and  3  and  for  all  the  calculations  of  figs*  8  to  10*  The 
plane-wave  values  were  also  determined  using  the  semi-empirical  methods 
of  French  ^).  The  factors  calculated  for  all  excited  states  of  ^B  up 
to  6*14  MeV  are  indicated  in  tables  5o  6,  and  7» 

The  first  point  to  note  is  the  similarity  between  the  plane-wave 
values,  and  many  of  the  EWBA  sets*  Although  the  absolute  spectroscopic 
factors  would  bear  no  relation,  the  relative  values  appear  to  be  compara¬ 
tively  stable,  and  all  the  sets  calculated  show  the  same  general  trends 
with  excitation  energy* 

The  absolute  EWBA  cross-sections  at  the  peak  angles  for  options  2 
to  5  where  the  same  potential  set  was  employed,  are  listed  in  table  8* 
Using  a  cut-off  of  3  fta  in  the  zero-range  approximation  would  appear  to 
reduce  the  absolute  cross-sections  by  a  factor  of  about  2,5  from  the 
finite-range  non-local  values*  However  the  spectroscopic  factors  for 
these  two  options  are  very  similar,  differing  by  no  more  than  15#  in 
most  cases* 

Also  in  all  cases,  the  inclusion  of  non-local  effects  tends  to 
slightly  increase  the  cross-section  over  the  values  determined  with  the 
finite-range  approximation  alone*  This  increase  occurs  despite  the 
fact  that  contributions  from  the  nuclear  interior  are  damped,  and  as 
already  mentioned,  is  due  to  the  non-local  enhancement  of  the  bound 
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TABLE  8.  ABSOLUTE  EWBA  CROSS-SECTIONS  (mb/sr)  at  0^.* 


State 

i 

ZRCO 

FRCO 

FR 

FRNL 

Gnd 

1 

7.45 

12.30 

17.60 

17.71 

0.72 

1 

3.25 

5.40 

7.80 

8.60 

1.74 

1 

1.48 

2.45 

3.30 

3.50 

2.1 5 

1 

5.05 

8.36 

11.00 

11.20 

3.58 

1 

12.00 

19.70 

24.30 

25.0 

4.77 

1 

22.40 

36.8 

43.00 

^4.50 

4.77 

2 

4.70 

6.30 

9.60 

9.70 

5.H 

0 

10.60 

13.70 

13.90 

15.60 

5.16 

1 

18.40 

30.10 

34.70 

35.70 

5.93 

1 

23.40 

38.20 

42.3 

43. 80 

6.14 

2 

6.77 

11.0 

11.2 

12.60 

♦Calculations  performed  using  only  the  set  of  parameters  in 
fig.  8. 
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The  average  EWBA  values  with  errors  are  compared  in  table  7  with 
1*  plane-wave  estimates 

2*  previous  measurements  of  Bus sine  and  Smith 
3*  theoretical  values  calculated  by  Kurath  ^). 

The  errors  were  estimated  from  the  spread  of  the  calculated  values*  An 
additional  10$  error  is  introduced  from  uncertainties  in  the  relative 
yields  to  the  various  states* 

The  agreement  with  theory  is  poor*  especially  for  the  0*72  MeV 
state,  and  in  general  the  average  values  do  not  closely  agree  with  the 
previous  measurements*  For  some  states  e»g*  0*72,  2*15  there  is  very 
close  agreement  between  the  averaged  and  the  plane-wave  factors,  while 
other  pairs  differ  by  almost  a  factor  of  2* 

4*4  Conclusions 

The  search  for  the  four  weak  states  of  has  resulted  in  nega¬ 
tive  evidence  for  all  of  them*  When  considered  in  the  light  of  previ¬ 
ous  measurements  using  this  and  other  reactions,  it  is  felt  that  the 
levels  at  5»58,  6*40,  6*77  MeV  are  spurious*  Although  the  5*18 
MeV  level  was  not  detected  in  this  reaction,  published  evidence  would 
appear  to  support  its  existence*  The  presence  of  a  level  at  this  exci- 

tation  neatly  explains  the  decay  of  the  7 <>56  MeV  level  ©f  B  )• 

10  10  ^ 

Also,  one  cannot  discard  the  results  from  the  B(d,d’)  B*  reaction  )* 
The  fact  that  no  significant  trace  has  been  found  would  indicate  that 
the  strong  5 *11  -  5*16  MeV  doublet  almost  completely  masks  whatever 
weak  excitation  is  present* 
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Q  ][Q 

From  the  IWBA  study  of  the  7Be(d,n)  B  reaction  at  5»5  MeV  several 
points  emerge *  The  first  point  is  that  in  any  stripping  reaction  it  is 
of  crucial  importance  to  determine  the  optical  parameters  to  be  later 
used  in  the  IWBA  analysis*  These  parameters  should  be  preferably  deter- 
mined  using  the  target  to  be  employed  in  the  strapping  reaction  study* 
(In  this  way,  isotopic  or  other  differences  should  not  affect  the 
parameters*) 

This  foreknowledge  will  avoid  the  somewhat  arbitrary  search  pro- 

ceedure  carried  out  in  this  report,*  In  previous  cases  where  parameters 

were  first  measured  (e*g*  ref*  JlJ)9  a  reasonable  degree  of  success  was 

achieved*  Even  for  this  classic  case  however,  several  parameter  sets 

yielded  impeccable  fits  to  the  scattering  data*  Unfortunately  these 

sets  resulted  in  differing  absolute  IWBA  cross-sections  and  hence 

differing  spectroscopic  factors*  This  was  also  the  case  in  the  calcu- 

6l 

lations  of  W*R*  Smith  )  who  concluded  that  spectroscopic  factors  are 
not  accurately  determined  using  the  present  IWBA  .formalism*  Despite 
these  difficulties,  prior  knowledge  of  the  parameters  is  by  far  the 
best  approach  to  IWBA  calculations* 

The  remarks  of  Schiffer  et  al  )  are  also  highly  relevant*  They 
point  out  that  with  present  day  developments  in  the  optical  model,  it 
is  not  sufficient  to  measure  optical  parameters  for  one  nucleus  at  one 
bombarding  energy  <,  Rather  one  should  measure  energy-  and  mass-averaged 
parameters  over  a  range  of  intermediate  energies,  with  masses  in  the 
region  of  the  nucleus  under  study* 
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The  second  point  of  interest  concerns  the  sensitivity  of  the  cal¬ 
culated  cross-sections  (and  hence  spectroscopic  factors)  to  optical 
parameter  uncertainties  and  the  precise  mod©  of  calculation.  The 
cross-sections  listed  in  table  8. for  a  single  parameter  set  illustrate 
how  drastically  cross-sections  can  change  with  the  different  options. 

In  a  sense,,  there  are  too  many  prescriptions  for  doing  DWBA  calcula¬ 
tions.  With  regard  to  this  point  Pinkston  and  Satchler  have  indi¬ 
cated  how  critically  the  calculations  depend  on  the  wave  function  used 
for  the  captured  particle.  In  carrying  out  a  search  for  DWBA  fits  to 
stripping  data  as  in  this  report ,  it  would  appear  best  to  confine  the 
calculations  to  the  FRNL  option.  This  represents  the  most  sophisti¬ 
cated  calculation  and  simulates  the  effects  of  using  a  radial  cut-off, 
by  dampening  contributions  from  the  nuclear  interior. 

The  efforts  to  reproduce  the  measured  angular  distributions  with 
one  set  of  optical  parameters  has  met  with  only  partial  success.  Gen¬ 
erally  the  main  stripping  peak  is  reproduced  for  the  excited  states  when 
using  one  set,  but  there  is  poor  agreement  with  the  cross-sections  at 
large  angles.  This  may  be  due  to  compound  nucleus  effects  which  will 
be  proportionately  more  important  at  these  angles. 

In  conclusion,  the  calculations  appear  to  show  that  at  this  energy, 
many  of  the  observed  transitions  in  the  ^Be(d,n)^^B  reaction  appear  to 
be  due  to  a  direct  reaction  mechanism  although  the  DWBA  theory  as  formu¬ 
lated  may  not  be  an  accurate  description. 
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APPENDIX  Ac  THE  TIME-SHIFT  CORRECTION 


The  origin  of  time-shift  is  easily  understood  if  one  considers 

36 

figo  A  )»  This  figure  shows  how  signal  timing  is  dependent  on 


pulse  amplitude  because  of  the  finite  discriminator  threshold*  The 
finite  bias  is  necessary  to  avert  triggering  on  random  noise  fluc¬ 
tuations. 

To  correct  for  this  time  shift  in  our  case,  the  amplitude  of  the 

detector 

signal  from  the  principal  neutron  is  recorded  by  the  computer 

simultaneously  with  the  neutron  time  signal*  Subsequently  a  channel 
shift  dependent  on  the  recorded  detector  signal  amplitude  is  applied 
to  the  neutron  count* 

The  appropriate  channel  shift  is  determined  prior  to  the  experi¬ 
ment  by  analysis  of  a  two-dimensional  4096  channel  record  of  the 
mono-energetic  neutron  group  from  a  reaction  such  as  T(p»n)^Hep  or 
D(d,n)^He.  One  dimension  corresponds  to  the  detector  signal  amplitude, 
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which  is  divided  into  thirty- two  groups  or  bins  of  pulse  height 9  the 
second  dimension  to  the  ordinary  time  or  energy  spectrum  of  the  neutrons. 
Thus  a  display  of  the  kQ9 6  data  channels  gives  thirty-two  time 
spectra  of  128  channels  each0  with  the  neutron  peak  successively  time 
-shifted  from  bin  to  bin0  CORRsa  subroutine  of  GPKS#  can  then  compute 
the  channel  shift  appropiate  to  each  bin  by  determining  the  centroid 
of  each  neutron  peak  relative  to  a  standard  peak0  Zero  channel  shift 
refers  to  the  bin  of  maximum  pulse  amplitudes. 

The  table  of  channel  shifts  is  then  stored  in  memory  and 
during  a  time-of- flight  experiment t  the  appropiate  channel  shift  is 
applied  to  each  time-event  recorded. 
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APPENDIX  B 


This  appendix  contains  relative  angular  distributions  for  the 
^Be(d,n)^B  and  ^0(dtn)^F  reactions  at  deuteron  bombarding  energy 
Ed  =  5*54  MeV.  There  are  fourteen  neutron  distributions  for  the 
^Be(d,n)^B  reaction  corresponding  to  ^B  excitations  up  to  7*00 
MeV.  The  two  distributions  for  the  “^0(d„n)^F  reaction  correspond 
to  excitation  of  the  ground  and  0.5  MeV  states  of  ^F. 
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ANGULAR  DISTRIBUTION  OF  THE  0.00  MeV  STATE  OF 
10 B  FROM  9Be(d„n)10B  AUG.  I8/65 


Angle  (C  of  M) 

Corrected  Area 

Error 

0.000 

0.00826 

0.00079 

5*539 

0.01121 

O.OOO85 

11.074 

0.01520 

0.00120 

16.601 

0. 0197? 

0.00161 

22.116 

0.02930 

0.00181 

27.615 

0.02778 

0.00192 

33*094 

0.02383 

0.00160 

38.550 

0.02068 

0.00148 

49*377 

0.01186 

0.00081 

65.362 

0.00739 

0.00059 

80.981 

O.OO65O 

0. 00049 

106.096 

0.00533 

0.00039 

115*815 

0. OO769 

0.00057 

125*357 

0.00956 

0.00062 

139*371 

0.00921 

0.00060 

153-089 

0.00891 

0.00065 

ANGULAR 

10 

DISTRIBUTION  OF  THE  0.717  MeV 

'B  FROM  9Be(d„n)10B  AUG.  I8/65 

STATE  OF 

Angle  (C  of  M) 

Corrected  Area 

Error 

0.000 

0.01229 

0.00097 

5*562 

0.01566 

0.00103 

11.120 

0.02026 

0.00141 

16.669 

0.02533 

0.00185 

22.206 

0.03655 

0.00207 

27*726 

0.03809 

0.00233 

33*225 

0.03008 

0.00183 

38.700 

0.02544 

0.00165 

49.563 

0.01111 

0.00075 

65.590 

0.00413 

0.00042 

81.236 

0.00389 

0.00038 

IO6.356 

0.00360 

0.00031 

H6.O63 

0.00373 

0.00039 

125*585 

O.OO363 

0.00037 

139*557 

0. 00354 

0.00037 

153*220 

0.00532 

0.00053 

0  STAta  •  iM  0O<0  3HT  50  i  31  Lf.TBIfl  fiAJl  Ak 
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ANGULAR  DISTRIBUTION  OF  THE  1.74  MeV  STATE  OF 
10 B  FROM  9Be(d»n)10B  AUG.  I8/65 


Angle  (C  of  M) 

Corrected  Area 

Error 

QoOQO 

0.00464 

0.00049 

5.600 

0.00528 

0.00052 

11.196 

0.00581 

0.00050 

16.783 

0,00589 

0.00058 

22.356 

0„008?8 

0.00068 

27.912 

0.00729 

0,00059 

33.446 

Oo 00564 

0.00055 

38.953 

0.00408 

0.00040 

49 0875 

0.00174 

0.00030 

65.973 

0.00114 

0.00023 

81.664 

0.00119 

0.00023 

106.792 

0.00052 

0.00017 

116.479 

0.00068 

0.00017 

125.968 

0.00084 

0.00021 

139.869 

0.00114 

0.00030 

153.441 

0.00086 

0.00027 

ANGULAR 

10 

DISTRIBUTION  OF  THE  2.15  MeV  STATE  OF 

B  FROM  ^B©(d„n)10B  AUG.  I8/65 

Angle  (C  of  M) 

Corrected  Area 

Error 

0.000 

0.00408 

0.00048 

5.618 

0.00474 

0.00048 

11.232 

0.00486 

0.00044 

16.836 

0.00514 

0.00054 

22.426 

0.00878 

0.00067 

27,998 

0.0088? 

0.00063 

33.548 

0.00746 

0.00061 

39.071 

0.00514 

0.00043 

50.020 

0.00444’ 

0.00041 

66.151 

0.00357 

0.00033 

81.862 

0.00319 

0.00031 

106.995 

0.00251 

0.00025 

116.672 

0.00360 

0.00028 

126.146 

0.00345 

0.00030 

140.014 

0.00455 

0.00047 

153.543 

0.00388 

0.00037 
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ANGULAR  DISTRIBUTION  OF  THE  3*58  MeV  STATE  OF 
10B  FROM  9B©(d!,n)1°B  AUG.  18/65 


Angle  (C  of  M) 

Corrected  Area 

Error 

OoOOO 

0*01304 

0*00080 

5*696 

0*01361 

0*00093 

11*386 

0*01276 

0.00093 

17.066 

0*0100? 

0.00086 

22*730 

0 *  01137 

0.00073 

28*3 75 

0.01077 

0*00071 

33 - 993 

0.00835 

0.00047 

39.582 

0*00719 

0.00052 

50*651 

0.00660 

0.00046 

66.92 5 

0*00566 

0*00047 

82*72? 

0*00516 

0*00033 

X07.877 

0*00452 

0*00029 

117.513 

0.00520 

0*00040 

126*921 

0*00566 

0*00035 

140*645 

0*00439 

0*00027 

153.988 

0*00412 

0*00037 

ANGULAR 

10 

DISTRIBUTION  OF  THE  4.??  MeV  STATE  OF 

'b  FROM  9Be(d0n)x0B  AUG.  I8/65 

Angle  (C  of  M) 

Corrected  Area 

Error 

0.000 

0*00546 

0*00065 

5.789 

0.00640 

0*00058 

11*572 

0*00500 

0*00059 

17.3*0 

0*00543 

0*00061 

23*09 7 

0.00507 

0*00052 

28*828 

0*00608 

0*00044 

34*530 

0*00585 

0*00044 

40*198 

O.OO55I 

0*00040 

51.432 

O.00630 

0*00050 

67.860 

0*00596 

0*00040 

83.772 

0*0065? 

0.0004? 

108*943 

0*00580 

0.00046 

118*529 

0*00638 

0*00057 

127.855 

0*00752 

0*00043 

141*406 

0*0048? 

0*00032 

154*525 

0*00293 

0.00037 
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ANGULAR  DISTRIBUTION  OF  THE  5«H  M©V  STATE  OF 
10 B  FROM  9Be(d„n)10B  AUG.  18/65 


Angle  (C  ©f  M) 

Corrected  Area 

Error 

OoOOO 

0,16293 

0.00379 

5  0823 

0,14526 

0.00482 

11 0  641 

0,08374 

0,00464 

1?  0  446 

O.05399 

0.00298 

23*233 

0,03440 

0.00179 

28*995 

0.01597 

0.00105 

34,728 

0,00834 

0.00070 

40^426 

0.00769 

0.00066 

51,694 

0©00976 

O.QQO83 

68.2Q6 

0,01053 

0.00058 

84*159 

0.00791 

0.00052 

109,337 

0.00301 

O.OOO38 

118*905 

0,00422 

0.00039 

128*201 

0,00458 

0,00041 

141,688 

0.00000 

0,00000 

154,723 

0.00183 

0.00018 

ANGULAR 

10 

DISTRIBUTION  OF  THE  5*16  MeV  STATE  OF 

'b  FROM  9Be(d9n)10B  AUG.  I8/65 

Angle  (C  ©£  M) 

Corrected  Area 

Error 

0,000 

0,02044 

0.00200 

5,829 

0.02168 

0.00212 

22,652 

0.03524 

0.00362 

17,462 

0.02818 

0.00242 

23,254 

0,01622 

0.00139 

29,022 

0,01127 

0.00097 

34,760 

0,00751 

0.00067 

40,462 

o, 00656 

0.00064 

51,738 

o, 00433 

0.00073 

68,261 

0,00171 

0.00034 

84,220 

0,00171 

0.00034 

109.400 

0.00262 

0.00037 

118,965 

0.00198 

0.00034 

128.256 

0,00696 

0.00068 

141,733 

0.00001 

0.00001 

154.755 

0.00159 

0.00041 

000,8 
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ANGULAR  DISTRIBUTION  OF  THE  5»93  MeV  STATE  OF 
9Be(d,n)10B  AUG.  18/65  GAUSS-FITTED 


Angle  (C  of  M) 

Corrected  Area 

Error 

0.000 

0.04674 

0.00371 

5.930 

0.04342 

0.00362 

11.852 

0.03860 

0.00321 

17.762 

0.03860 

0.00306 

23.650 

0.03908 

0.00310 

29.512 

0.02870 

0.00184 

35.340 

0.01953 

0.00125 

41.129 

0.01397 

0.00146 

52.562 

0.00730 

0.00060 

69.274 

0.00462 

0.00045 

85.354 

0.00497 

0.00042 

110.556 

0.00256 

0.00039 

120.066 

o. 00279 

0.00035 

129.269 

o. 00357 

0.00044 

1^2.557 

0.00389 

0.00058 

155.335 

0.00348 

0.00056 

ANGULAR  DISTRIBUTION  OF  THE 

6.02 

MeV  STATE  OF 

ioB 

FROM  9Be(d9n)10B 

AUG. 

18/65 

Angle  (C  of  M) 

Corrected  Area 

Error 

0.000 

0.00760 

0.00085 

5.944 

0.00677 

0.00097 

11.881 

0.00479 

0.00066 

17.805 

0.00450 

0.00082 

23. 708 

0.00464 

0.00077 

29.583 

0.00435 

0.00088 

35.424 

0.0029 6 

O.OOO69 

41.225 

0.00164 

0.00082 

52.682 

0.00176 

0.00045 

69.421 

0.00286 

0.00044 

85.518 

0.00001 

0.00001 

110.723 

0.00163 

O.OOO36 

120.226 

0.00192 

0.00033 

129.416 

0.00200 

0. 00037 

142.676 

0.00142 

0.00043 

155.419 

0.00034 

0.00033 
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ANGULAR  DISTRIBUTION  OF  THE  6.14  MeV  STATE  OF 
10B  FROM  9Be(d,n)10B  AUG«  I8/65 


Angle  (C  of  M) 

Corrected  Area 

Error 

0.000 

0.0117 

0.0006 

5o965 

0.0112 

0.0006 

11.923 

0.0118 

0.0007 

17.867 

0.0132 

0.0007 

23.790 

0.0132 

0.0006 

29.685 

0.0132 

0.0008 

35.545 

0.0115 

0.0006 

41.364 

0.0124 

0.0012 

52.854 

0.0097 

0.0006 

85.755 

0.0045 

0.0004 

IIO.965 

0.0029 

0.0004 

120.456 

0.0026 

0.0004 

129.628 

0.0029 

0.0005 

142.848 

0.0020 

0.0005 

155.540 

0  0  00X2 

0.0004 

ANGULAR 

DISTRIBUTION  OF  THE 

6.57  MeV  STATE 

OF 

10B  FROM  ^Be(d9n)10B 

AUG.  18/65 

Angle  (C  of  M) 

Corrected  Area 

Error 

0.000 

0.0078 

0.0022 

6.050 

0.0058 

0.0023 

12.092 

0.0063 

0.0015 

18.118 

O.OO93 

0.0024 

24.122 

0.0127 

0.0011 

30.096 

O0OI56 

0.0009 

36.032 

0.0156 

0.0009 

41.924 

0.0168 

0.0009 

53.546 

0.0140 

0.0009 

70.485 

0.0071 

0.0011 

86.709 

O.OO3I 

0.0004 

111.939 

0.0052 

0.0017 

121.383 

0.0000 

0.0002 

130.480 

0.0025 

0.0008 

143.540 

0.0001 

0.0014 

156.027 

0.0000 

0.0006 

' 
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ANGULAR  DISTRIBUTION  OF  THE  6*88  MeV  STATE  OF 
10 B  FROM  9B©(d9n)10B  AUG.  18/65 


Angle  (C  of  M) 

Corrected  Area 

Error 

0.000 

0.0025 

0.0008 

6.1 2? 

0.0029 

0.0006 

12.24? 

0.0017 

0.0006 

18.350 

0.0014 

0.0006 

24.428 

0.0011 

0.0004 

30.474 

0.0012 

0.0003 

36.481 

0.0015 

0.0002 

42.439 

0.0001 

0.0001 

54.184 

0.0009 

0.0003 

71.271 

0.0005 

0.0002 

87.590 

0.0000 

0.0000 

112.838 

0.0000 

0.0000 

122.239 

0.0001 

0.0002 

131.266 

0.0000 

0.0000 

144.178 

0.0000 

0.0000 

ANGULAR  DISTRIBUTION  OF  THE  7.01 

MeV  STATE  OF 

10 B  FROM 

9Be(d„n)10B  AUG. 

18/65 

Angle  (C  of  M) 

Corrected  Area 

i 

Error 

0.000 

0.0049 

0.0009 

6.166 

0.0035 

0.0006 

12.323 

0.004? 

0.0007 

18.463 

0.0044 

0.0007 

24.578 

0.0023 

0.0004 

30.660 

0.0013 

0.0003 

36.701 

0.0009 

0.0002 

42.692 

0.0010 

0.0003 

54.498 

0.0000 

0.0000 

71.658 

0.0002 

0.0002 

88.024 

0.0000 

0.0000 

113.280 

0 . 0001 

0.0001 

122.661 

0.0000 

0.0000 

131.653 

0.0000 

0.0000 

144.492 

0.0000 

0.0000 
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ANGULAR  DISTRIBUTION  OF  THE  0.00  MeV  STATE  OF 
l7F  FROM  l60(d„n)17F  DEC,  8/65 


Angle  (C  of  M) 

Corrected  Area 

Error 

OoOOO 

0,1181 

0,0059 

5*530 

0,1259 

0,0067 

11.0 57 

0,1680 

0,0088 

I60 575 

0,2267 

0,0118 

22 . 082 

0,2686 

0,0140 

2 7-572 

0,3133 

0,0163 

33 0 044 

0 „ 3285 

0,0172 

38,492 

0,3058 

0,0160 

49 0 306 

0,22^2 

0,0119 

65,276 

0,1014 

0,0056 

80 0886 

0,0446 

0,0026 

96,094 

O0O35O 

0,0021 

IO60 001 

0 , 027 0 

0,001? 

115,725 

0,0522 

0,0031 

125,274 

0,0580 

0,0035 

139,304 

0,0537 

0,0033 

153 0 042 

0,0532 

0,0033 

5,530 

0,1260 

0,0067 

ANGULAR 

DISTRIBUTION  OF  THE  0,50  STATE  OF 

1?F  FROM  1^0(dttn)1^F  DEC,  8/65 

Angle  (C  of  M) 

Corrected  Area 

Error 

0,000 

1,2908 

0,0690 

5,576 

1.1861 

0,0634 

11,148 

0,9084 

0,0486 

16,711 

0,5937 

0,0319 

22,262 

0,3206 

0,0173 

27,795 

0,1534 

0,0084 

33,307 

0,0769 

0,0043 

38,794 

0,0549 

0,0031 

49,679 

0,0884 

0.0049 

65,734 

0,1010 

0,0058 

81,397 

0,0725 

0.0043 

96,623 

0,0372 

0,0023 

106,522 

0,0180 

0.0012 

116,222 

0,0276 

0,0018 

125,732 

0,0221 

0.0015 

139,677 

0,0184 

0.0013 

153,305 

0,0141 

0,0010 

5,576 

1,1314 

0.0606 

_  _ _ _  .  . 
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APPENDIX  C 

The  use  of  the  time-of-f light  system  to  determine  target  thickness. 

The  feasibility  of  estimating  target  thickness  by  means  of  the 
neutron  time-of-flight  spectrometer  has  been  investigated.  To  this 

-j  £  i  n 

end  thin-  and  thick-target  yields  of  neutrons  from  the  0(d,ni)  'F 

transition  were  studied  at  3® 50  MeV  deuteron  energy. 

The  thin  target  has  already  been  described  in  section  2.4,  and 

q 

consisted  of  a  thin  film  of  Be  deposited  on  a  Molybdenum  foil.  The 

^0  contamination  is  evident  in  the  ^Be(d,n)^B  spectrum  of  fig.  3® 

o  17 

At  0  ,  the  neutron  peak  corresponding  to  the  first  excited  state  of  ¥ 

(0.5  MeV  excitation)  is  not  contaminated  by  neutrons  corresponding  to 

states. 

The  thick  target  consisted  of  a  small  foil  of  Zirconium  canning 
material  used  in  one  of  the  Chalk  River  reactors.  The  surface  of  this 
material  had  been  oxidised  by  Oxygen-16  and  Oxygen  17 5  and  a  request 
had  been  made  for  some  estimate  of  the  depth  of  Oxygen  penetration 
into  the  Zirconium  surface. 

If  the  electronic  resolution  of  the  time-of-flight  spectrometer  is 

l6  17 

very  good,  then  the  shape  of  the  neutron  peak  for  the  0(d,n^)  F  reaction 

l6 

should  reflect  the  distribution  of  the  0  in  the  target.  The  fastest 
neutrons  will  come  from  the  surface  layers  of  the  target,  and  will  hence 
have  the  shortest  flight  times.  As  the  deuterons  penetrate  the  target 
they  lose  energy  by  ionisation  of  the  target  atoms  and  will  thus  give  rise 
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FIGC-2.  THIN-TARGET  YIELD  OF  YIELD  OF  NEUTRONS 

FROM  THE  l60(d,n,)l7F  TRANSITION  AT  3-5 
MEV  DEUTERON  ENERGY. 
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to  neutrons  of  progressively  lower  energies  as  the  depth  of  penetration 
increases.  The  net  effect  of  this  deuteron  energy  loss  is  to  broaden  the 
peak  in  the  neutron  spectrum. 

The  broadened  shape  is  a  function  of  both  the  electronic  resolution 
and  the  distribution  of  Oxygen-16  in  the  target.  Since  the  deuteron 
absorption  cross-sections  for  Oxygen  and  Zirconium  are  small  (<  1  barn) 
we  can  assume  that  a  uniform  number  of  deuterons  interacts  with  all  layers 
of  the  target.  (The  absorption  cross-sections  for  deuterons  of  this  energy 
were  estimated  using  Perey's  optical  model  code.) 

The  experiment  was  carried  out  using  the  two  targets  and  the  time- 
of- flight  system  described  in  Chapter  II.  An  average  current  of  approx¬ 
imately  lvA  of  MeV  deuterons  was  employed  and  the  principal  neutron 

detector  was  fixed  at  a  flight  path  of  6.3m  and  a  laboratory  angle  of  0°. 

Unfortunately  the  thick  target  was  seriously  affected  by  the  deuteron 
beam  and  only  one  or  two  satisfactory  thick  target  yields  were  collected. 
Thick-  and  thin-target  yields  are  illustrated  in  figs.  C-2,  and  C-l. 

The  neutron  energy  corresponding  to  this  transition  is  given  kinem¬ 
atically  as  1.229  MeV  where  the  Q-value  is  -2.131  MeV.  The  full-width  at 
half-maximum  (FWHM)  for  the  thin  target  is  2. 37  ns„  The  latter  peak  is  almost 
Gaussian  in  shape  but  has  a  'tail*  on  the  low  energy  side  of  the  peak. 

The  thick  target  peak  is  broadened  considerably  to  about  8ns  FWHM.  It 

l6 

is  evident  from  the  asymmetric  shape  that  the  distribution  of  0  in  the 
target  is  a  non-linear  function  of  the  depth.  If  the  "^0  had  been  uniform, 
a  roughly  rectangular  peak  shape  should  result. 
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If  the  delays  introduced  by  the  target  energy  losses  are  negli¬ 
gible  then  the  neutron  peak  represents  the  intrinsic  resolution  curve 
of  the  electronics,  P(t),  where  t  is  the  flight  time*  If  there  is  a 
measurable  energy  loss  due  to  target  thickness,  then  one  obtains  the 
delayed  curve  F(t)*  One  can  write 

F(t)  =/°°f(tJ)oP(t-t5)dtJ 

•=»  OO 

where  f(t8)  is  related  to  the  distribution  of  D0  atoms  in  the  target 
associated  with  flight  time  t" ,  and  this  function  is  the  object  of 
interest* 

If  we  can  assume  that  P(t)  is  given  by  the  thin  target  curve,  then 
we  can  take  some  functional  form  for  f(t),  and  numerically  integrate  to 
find  F(t)*  As  a  first  approximation,  it  was  decided  to  try  an  exponen¬ 
tial  distribution  for  the  target  atoms.  Thus  the  number  of  ^0  atoms 
in  (x,x+dx)  is  const* exp (-o(x)dx  where  x  is  the  depth  coordinate  measured 
normal  to  the  target  surface*  For  the  reaction  under  study,  the  depth  x, 
and  the  flight  time  are  related  kinematically  as 

x  =  2*27*10^  l/to  -  l/t2!  cm 

|dEc(/dx|  '  / 

where  t0  is  the  flight  time  in  ns  for  a  1*229  MeV  neutron  over  a  6. 3m 
flight  path  (=410* 5ns) ? 

t  is  the  flight  time  in  ns  for  a  neutron  coming  from  depth  x  in  the  target; 
|dEd/dx|  is  the  tabulated  target  energy  loss  in  units  of  MeV/cm*  (This 
quantity  has  been  assumed  constant  in  deriving  the  above  formula*) 
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The  tables  of  Coon  (ref*  C)  gave 

|dEd/dx|=  101  keVorag^.cm 

Thus 

exp(-<£c)dx  =  exp  (-2.21.o(.10~3.(l  -  l/(l  +  t/410.5)2))dt 

(410.5  +  t)3 

For  convenience  it  was  decided  to  take  a  Gaussian  form  for  P(t)  with 
a  standard  deviation  estimated  from  curve  C-l*  This  leads  to 

roo  p 

F(t)  =  /  exp(-0.466(t-t'r  -  2.21,°(.10~3.(1  -1/(1  +  t'/410.5)  ))dt' 

Jo  (410.5  +  t')3 

The  distribution  function  fCt”)  =  0  for  t*<  0„  and  this  condition  fixes 
the  integration  limits* 

The  above  integral  for  F(t)  was  computed  numerically  for  several 
values  of  o( ,  until  a  reasonable  fit  to  the  asymmetric  experimental  shape 
was  obtained*  The  results  of  one  such  calculation  for  c<=  9300cm”^  are 
indicated  in  fig*  C-2*  Acceptable  agreement  with  the  shape  could  be 
obtained  for  values  of  o<_  within  ±  of  the  above  „  at  least  for  about  60 
channels  from  the  peak  channel*  Thereafter  the  theoretical  shape  under¬ 
estimated  the  experimental  counts*  This  discrepancy  may  be  due  to  any 
one  of  a  number  of  effects  including  contributions  from  the  tail  of  the 
experimental  P(t)  or  breakdown  of  the  assumption  of  an  exponential  dis¬ 
tribution  of  target  atoms*  However  the  agreement  over  most  of  the  curve 
tends  to  support  the  latter  assumption,,  which  can  also  be  justified  on 
physical  grounds  if  the  target  was  formed  by  a  diffusion  of  ^0  atoms  into 
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a  pure  Zirconium  surface.  With  the  estimated  value  of  cxC (9300cm”'*') 
the  number  of  target  atoms  drops  to  half  the  number  at  the  surface  at 
a  depth  of  about  7500  angstroms. 

In  conclusion  this  rough  analysis  appears  to  show  that  it  is 
possible  to  obtain  some  idea  of  the  distribution  of  target  atoms  from 
the  shape  of  the  associated  neutron  peak.  As  a  method  of  estimating 
target  thickness  it  holds  promise  although  more  experimentation  and 
analysis  is  necessary  to  place  the  technique  on  a  firmer  basis. 
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